DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 

Wright-Patterson  Air  Force  Base,  Ohio 

89  12  14  0^4 


STATEMiNT  A 
»d  for  puEIic  r«!*cu«; 


AFIT/GSO/ENS/89D-10 

IMAGE  NAVIGATION  OF  TIROS-N 
WEATHER  SATELLITE  DATA 

THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Space  Operations 


Charles  H.  Larcomb,  B.A. 

Captain,  USAF 

December  1989 

Approved  for  public  release;  distribution  unlimited 


Preface 


Hobbyists  have  been  receiving  and  displaying  weather 
satellite  images  for  some  time.  With  relatively  sim.ple  and 
inexpensive  ground  stations,  enthusiasts  can  receive  Automa¬ 
tic  Picture  Transmissions  (APT)  from  both  the  U.S.  TIROS-N 
satellites  and  the  Soviet  Meteor  satellites.  The  TIROS-N 
polar-orbiting  satellites  continuously  transmit  weather  data 
as  they  orbit  the  earth  every  102  minutes. 

The  main  obstacle  to  using  these  images  is  a  lack  of 
geographic  references.  Geographic  boundaries  are  often 
obscured  by  clouds  or  lack  of  contrast.  This  work  is  aimed 
at  removing  this  obstacle,  to  greatly  enhance  the  usefull- 
ness  of  polar-orbiting  weather  satellite  images. 

I  am  indebted  to  my  thesis  advisor,  Dr.  Thomas  S. 

Kelso,  for  his  guidance  and  inspiration  through  no  small 
expense  of  time  and  personal  resources.  Ms.  Mona  Smith  of 
NOAA/NESDIS  was  very  helpful  in  obtaining  the  NOAA  publica¬ 
tions  used  in  this  work.  Finally,  I  thank  my  wife,  Cirela, 
for  her  patience  and  support. 

Charles  H.  Larcom.b 
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Abstract 

The  purpose  of  this  study  was  to  develop  an  algorithm 
to  perform  spatial  registration  of  (assign  geographic  coor¬ 
dinates  to)  TIROS-N  weather  satellite  data  using  a  personal 
computer.  Specific  objectives  were  (1)  to  gather  currently 
available  information  into  a  unified  and  understandable 
form,  (2)  to  develop  the  background  and  related  theory 
necessary  for  the  registration  of  both  High  Resolution 
Picture  Transmission  (HRPT)  and  Automatic  t'lcture  Transmis¬ 
sion  (APT)  data,  and  (3)  to  present  a  specific  registration 
algorithm  for  TIROS-N  weather  satellite  data. 

The  approach  chosen  is  to  use  a  general  elliptical 
orbit  model  in  conjunction  with  accurate  satellite  element 
sets  (orbital  parameters)  to  calculate  accurate  satellite 
positions  at  a  given  time.  Knowing  the  earth  movement  and 
sensor  operation  as  functions  of  time,  the  geographic  loca¬ 
tions  of  each  pixel  (picture  element)  are  determined  through 
geometric  relations. 

The  algorithm  presented  converts  pixel  coordinates  to, 
geographic  coordinates  (direct  referencing)  and  vice  versa 
(inverse  referencing).  These  algorithms  assume  an  ellipti¬ 
cal  orbit  (with  perturbations)  and  an  oblate  spheroid  earth. 


xiii 


The  most  critical  factor  is  found  to  be  timing.  A 
timing  error  of  one  second  in  computing  satellite  position 
produces  an  earth  location  error  of  about  6.5  kilometers  at 
the  satellite  subpoint.  A  timing  correction  method  is 
outlined  to  visually  improve  the  image  registration. 
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SPATIAL  REGISTRATION  OF  TIROS-N 
WEATHER  SATELLITE  DATA 

I .  Introduction 

Background 

The  TIROS-N  satellite  series  is  the  set  of  polar  orbit¬ 
ing  satellites  operated  by  the  National  Oceanographic  and 
Atmospheric  Administration  (NOAA)  that  provides  advanced  en¬ 
vironmental  (weather)  data.  The  active  satellites  con¬ 
tinuously  transmit  both  High  Resolution  Picture  Transmission 
(HRPT)  data  (1.1  km  resolution)  and  Automatic  Picture  Trans¬ 
mission  (APT)  data  (4  km  resolution)  as  they  orbit  the  earth 
at  an  altitude  of  about  850  kilometers  [1:2-1, 3-5]. 

The  images  do  not  contain  geographic  grids  or  map 
outlines  as  provided  by  the  more  sophisticated  high-altitude 
Geosynchronous  Operational  Environmental  Satellites  (GOES) . 
These  enhancements  must  be  added  manually  after  image  recep¬ 
tion.  This  spatial  registration  (process  of  assigning  geo¬ 
graphic  coordinates  to  the  satellite  data)  is  essential  to 
provide  accurate  and  detailed  weather  information. 

The  HRPT  signal  has  a  high  data  rate  and  frequency.  Its 
reception  requires  a  fairly  sophisticated  ground  station 
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using  mainframe  computers  which  also  allow  partial  automa¬ 
tion  of  the  registration  process.  Current  HRPT  ground 
stations  are  expensive,  costing  $150,000  or  mere  [17:2-2]. 

The  APT  signal  has  a  low  data  rate  and  frequency  per¬ 
mitting  reception  with  relatively  simple  and  inexpensive 
receivers  [1:55].  A  complete  ground  station  can  be  con¬ 
structed  for  less  than  $500  [17:4-1].  These,  however,  gen¬ 
erally  provide  film  or  paper  copies  which  must  be  regist¬ 
ered  by  hand  using  tedious  and  inaccurate  methods. 

With  the  growing  popularity  of  personal  computers, 
computer  display  systems  are  becoming  a  common  method  of 
displaying  weather  satellite  images  [17:9-4].  For  a  modest 
cost,  an  analog-to-digital  converter  allows  PCs  to  process 
and  store  APT  satellite  image  data.  Commercially  available 
software  is  expensive,  though,  and  spatial  registration  is 
still  a  problem. 

Due  to  these  difficulties,  the  usefulness  of  TIROS-N 
weather  data  has  been  largely  unexploited.  Most  weather 
agencies  depend  heavily  on  data  from  a  few  costly  GOES 
satellites  of  which  only  one  is  currently  operating. 

TIROS-N  data  with  spatial  registration  could  provide  an 
essential  and  capable  backup  capability,  as  well  as  coverage 
of  polar  regions.  Simple,  mobile  APT  receivers  could  also 
be  easily  transported  to  any  point  on  the  globe. 


Problem  Statement 


The  purpose  of  this  research  is  to  develop  an  algo¬ 
rithm  to  register  TIROS-N  weather  satellite  data  on  a  per¬ 
sonal  computer.  After  the  analog  (continuous)  APT  data  from 
the  satellite  is  reconverted  to  digital  form,  the  algorithm 
will  assign  geographic  coordinates  to  each  pixel  (individual 
data  element)  or  find  the  pixel  corresponding  to  a  given 
geographic  location.  This  registration  allows  the  user  to 
add  latitude/longitude  grids  and  map  overlays,  display, 
process,  analyze,  and  store  the  satellite  images  on  a  per¬ 
sonal  computer. 

Objectives 

This  research  was  directed  to  meet  the  following 
objectives: 

1.  Gather  currently  available  information  necessary 
for  HRPT  and  APT  data  registration  into  an  acces¬ 
sible,  unified,  and  understandable  form 

2.  Develop  the  background  and  related  theory  for  the 
spatial  registration  and  use  of  both  HRPT  and  APT 
data ;  and 

3.  Present  a  specific  registration  algorithm  for 
TIROS-N  weather  satellite  data 
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II .  Literature  Review 

Introduction 

As  one  author  recently  observed  [6:898],  "there  are  few 
published  articles  on  the  image  referencing  and  geometric 
correction  of  AVHRR  data."  A  current  literature  search 
confirmed  rhis  deficit  and  uncovered  no  articles  on  the 
registration  of  APT  data. 

Nevertheless,  the  basic  principles  that  apply  to  spa¬ 
tial  registration  are  well  known.  The  methods  used  for  the 
higher  resolution  HRPT  data  should  also  apply  in  large  part 
to  APT  data. 

The  following  paragraphs  review  literature  on  the 
topics  pertaining  to  this  investigation.  These  topics  are 
spacecraft  and  sensor  characteristics,  satellite  orbits, 
geometric  corrections,  HRPT  registration  methods,  and  AVHRR 
data  applications. 

Spacecraft  and  Sensor  Characteristics 

The  TIROS-N  series  currently  consists  of  two  operation¬ 
al  satellites  in  sun-synchronous  polar  orbits  [8:1-7]. 
Detailed  spacecraft  and  sensor  characteristics  are  presented 
in  the  NOAA  User's  Guides  [1;  8].  These  include  orbital 
characteristics,  specific  parameters  on  spacecraft  instru- 


merits,  and  the  formats  of  archive  data  (stored  data  avail¬ 
able  from  NOAA  upon  request) . 

The  heart  of  the  TIROS-N  series  spacecraft  sensors  is 
the  Advanced  Very  High  Resolution  Radiometer  (AVHRR) .  This 
instrument  provides  real-time  visible  and  infrared  data 
worldwide  in  two  formats:  Automatic  Picture  Transmission 
(APT)  with  4  km  resolution,  and  High  Resolution  Picture 
Transmission  (HRPT)  with  1.1  km  resolution  [1:2-1]. 

NOAA  Technical  Memorandum  NESS  95  [15]  provides  detail¬ 
ed  information  on  the  operation  of  the  AVHRR  and  other 
sensors  and  satellite  systems.  The  format  and  parameters 
for  both  HRPT  and  APT  data  transmissions  is  also  included. 

TM  9  5  was  the  major  source  of  background  infoirmation  presen¬ 
ted  in  this  work. 

NOAA  Technical  Memorandum  NESS  107  -  Rev  1  [10]  provi¬ 
des  information  useful  for  the  calibration  (assigning  more 
accurate  absolute  measurement  values)  of  TIROS-N  sensor 
data.  Also,  Appendix  C  outlines  a  generic  registration 
algorithm  for  polar-orbiting  satellites. 

Satellite  Orbits 

Ideally,  a  satellite  in  a  circular  orbit  about  the 
earth  follows  a  perfect  circle  determined  by  the  laws  of 
orbital  mechanics  for  two-body  motion.  In  reality,  however, 
there  are  deviations  due  to  a  slight  eccentricity  of  the 
satellite  orbit,  and  differential  gravity  effects  caused  by 


the  sun  and  moon,  and  an  earth  that  bulges  at  the  equator. 
Measuring  the  exact  time  is  also  critical  for  accurate 
position  determination. 

Two  approaches  have  been  used  to  tackle  these  problems. 
One  [6:898;  11:13]  assumes  an  ideal  circular  orbit.  Cor¬ 
rections  are  then  made  based  on  known  deviations  such  as 
orbit  eccentricity,  or  using  known  image  features. 

An  alternate  method  is  to  use  the  actual  orbital  param¬ 
eters  in  conjunction  with  precise  time  marks  included  in  the 
data  transmissions  [18:1258].  This  method  yields  higher 
accuracy  than  the  former,  but  requires  timely  and  accurate 
measurements  of  the  satellite  orbital  parameters.  This  is 
the  approach  used  in  this  work. 

Geometric  Corrections 

Once  the  spacecraft  position,  sensor  attitude,  and 
earth  position  are  known  at  a  given  time,  the  relation 
between  an  observed  image  point  and  its  corresponding  point 
on  the  earth's  surface  is  a  matter  of  geometry.  This  geome¬ 
try  includes  distortions,  however,  which  must  be  corrected 
to  accurately  display  satellite  images.  The  distortions  due 
to  a  rotating  earth  and  moving  satellite  during  an  image 
scan  are  relatively  straightforward.  The  view  of  the 
earth's  surface  from  an  angle  is  more  complicated,  since  the 
earth  is  not  a  perfect  sphere.  Because  the  satellite  is  in 
low  orbit  (about  850  km)  with  a  wide  field  of  view  (+  55.4 
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degrees) ,  this  skewed  view  of  the  earth  causes  the  greatest 
distortions.  [18:1257]. 

There  are  two  conventional  distortion  correction  metho¬ 
ds.  One  is  to  use  a  spherical  model  of  the  earth  and  make 
corrections  for  local  deviations.  This  method  yields  errors 
of  several  pixels  and  runs  slowly.  An  alternate  method  uses 
known  image  reference  points  called  Ground  Control  Points 
(GCPs) .  While  very  precise,  this  method  depends  on  having  a 
number  of  clear  and  identifiable  GCPs  in  the  image 
[18:1257].  This  method  may  be  better  for  registration 
verification. 

Sun  [18:1257]  proposes  an  improved  method  based  on  an 
elliptical  model.  Computations  with  this  method  are  four 
times  faster  than  the  spherical  model  with  corrections,  and 
have  less  than  three  pixel  errors  [18:1257]. 

The  approach  taken  in  this  work  is  to  provide  an  accur¬ 
ate  image  registration  based  on  a  general  elliptical  satel¬ 
lite  orbit  model  and  an  oblate  spheroid  earth  model.  By 
accurately  determining  the  geographic  location  of  each 
pixel,  images  can  be  displayed  without  distortion  in  any 
format  desired. 

HRPT  Registration  Methods 

A  few  authors  propose  specific  methods  for  registering 
AVHRR  HRPT  data.  Popta  [11]  uses  a  method  to  superimpose 
geographic  data  (latitude/longitude  crosses  and  map  out- 


7 


lines)  on  HRPT  images  in  real-time.  The  model  used  assumes 
a  circular  orbit  about  a  spherical  earth,  with  corrections 
to  account  for  a  rotating,  oblate  earth,  and  variable  satel¬ 
lite  speed,  eccentricity,  equator  crossing  time,  and  nodal 
period.  Errors  of  10  km  at  subpoint,  and  up  to  21  km  at 
edges  are  noted,  due  to  "low  accuracy"  orbital  data  [11:20- 
21]  . 

This  method  is  near  real-time  (5  minutes  run  time 
before  data  receipt),  but  has  several  disadvantages.  Since 
the  orbital  track  is  different  for  each  pass,  images  cannot 
be  compared  directly  in  sec[uence  for  time-series  analysis. 
This  method  does  not  use  the  accurate  time  marks  present  in 
the  data.  Also,  geographic  coordinates  are  only  assigned  to 
image  locations  where  the  grid  lies,  not  to  all  image  loca¬ 
tions  [11:1,3,9] . 

Tozawa  presents  another  registration  approach.  Here, 
AVHRR  data  is  mapped  to  a  surface  coordinate  system  using  a 
mapping  function.  This  function,  solved  by  approximation, 
is  based  on  current  satellite  orbit  parameters  provided  by 
NOAA.  GCPs  are  not  used  in  this  method  because  the  mapping 
function  is  very  complex  over  the  3000  km  image  swath. 

Also,  few  GCPs  can  be  found  in  ocean  areas. 

The  author  claims  the  procedure  is  fast  and  accurate. 
However,  he  uses  one  GCP  to  adjust  the  data  presented, 
citing  inaccurate  orbital  data  [19:46,51]. 
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Ho  [6]  provides  a  better  method.  A  corrected,  spheri¬ 
cal  model  and  a  spherical  satellite  orbit  are  used.  Nodal 
information  and  the  time  of  first  scan  line  are  included. 

The  apparent  satellite  height  and  inclination  are  adjusted 
using  one  GCP.  Finally,  the  resulting  registration  is 
verified  using  11-25  GCPs  for  reference.  Average  errors  are 
8  pixels  without  the  GCP,  and  3  pixels  with  it  [6:898-900]. 

The  error  rate  with  one  GCP  is  excellent.  A  disad¬ 
vantage  for  real-time  processing  is  that  the  operator  has  to 
interrupt  the  registration  to  locate  a  GCP.  While  the 
author  asserts  that  GCP  location  can  be  automated,  no  method 
is  suggested. 

An  excellent  registration  method  is  presented  by  Sun 
[18]  using  a  two-step  approach.  First,  an  ellipsoidal  model 
of  the  earth  provides  corrections  for  the  entire  image.  The 
time  code  imbedded  in  each  scan  line  is  used  to  determine 
the  satellite  position.  Second,  a  template  matching  using 
6-10  GCPs  corrects  residual  distortions  in  a  local  area  of 
interest . 

The  first  step  reduces  errors  to  0-3  pixels  and  works 
four  times  faster  than  a  spherical  model  with  corrections. 

If  greater  accuracy  is  desired,  the  second  step  removes 
virtually  all  geometric  distortions  [18:1257]. 

Sun  also  employs  an  interesting  time  reduction  techni¬ 
que  to  process  over  2048x4500  pixels  present  in  each  image. 
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4x17  pixel  blocks  are  processed  by  performing  calculations 
for  the  center  point  and  using  bilinear  interpolation  for 
the  remaining  points  in  the  block.  While  not  an  objective 
of  this  thesis,  this  method  could  be  applied  to  reduce 
computer  processing  time. 

Emery  [3]  provides  a  good  summary  and  review  of  current 
AVHRR  registration  techniques.  A  sample  application  using 
HRPT  data  compares  their  relative  performance. 

AVHRR  Data  Applications 

Besides  the  obvious  use  of  AVHRR  data  for  routine  weather 
observation  and  forecasts,  there  are  other  useful  applica¬ 
tions.  Compared  to  LANDSAT  images  AVHRR  data  has  signi¬ 
ficant  advantages,  despite  the  reduced  resolution  [5;  9]. 
AVHRR  images  are  available  2-4  times  each  day,  compared  to 
once  every  18  days  (at  best)  for  LANDSAT.  The  field  of  view 
is  much  wider  for  AVHRR,  providing  greater  coverage  [5: GI¬ 
GS].  These  advantages  proved  to  be  an  attractive  tool  for 
agricultural  monitoring  over  France  [7:1041]  and  obser¬ 
vation  of  sea  surface  temperatures  near  Japan  [9:1259]. 

While  APT  data  has  lower  resolution  than  HRPT  data,  it 
could  prove  to  be  very  useful  because  of  its  widespread 
availability.  Unlike  HRPT,  receipt  of  APT  is  not  limited  to 
a  few  expensive  ground  receiving  stations. 

The  approach  of  this  research  is  to  use  accurate  posi¬ 
tion  calculations  with  precise  time  code  information  for 
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image  registration.  This  should  allow  fully  automated 
operation  and  still  provide  an  accurate  registration. 

Summary 

The  basic  principles  involved  in  image  registration  are  well 
known.  While  no  application  has  been  made  to  APT  data,  a 
few  methods  have  been  developed  for  HRPT  data.  Due  to  the 
data  similarities,  these  provide  an  excellent  approach  for 
APT  data. 

The  problem  presented  involves  a  trade-off  between 
processing  time  and  registration  accuracy.  Sun  presents  an 
excellent  two-step  method  which  seems  to  optimize  this 
trade-off,  at  least  for  HRPT  data. 

The  two  basic  registration  methods  are  to  use  an  accur¬ 
ate  orbital  prediction  model,  or  to  manually  identify  known 
GCP  earth  locations.  The  approach  used  here  is  a  general 
elliptical  orbit  model  (with  perturbations) ,  an  oblate 
spheroid  earth  model,  and  accurate  time  code  information. 

AVHRR  data  has  important  applications  in  addition  to 
weather  monitoring.  The  development  of  a  registration 
algorithm  for  both  HRPT  and  APT  data  will  make  TIROS-N 
weather  data  even  more  useful  and  accessible. 
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Introduction 


The  earth  data  viewed  by  the  satellite  consists  of  a 
series  of  individual  picture  elements  (pixels) .  Like  tiles 
on  a  mosaic,  these  form  an  image  of  the  earth  scene  as  the 
satellite  sensor  scans  the  area. 

The  location  on  the  earth's  surface  viewed  in  each 
pixel  is  determined  by  (1)  the  satellite  position  and  orien¬ 
tation  in  space  (2)  the  sensor  operation,  and  (3)  the 
earth's  orientation  in  space.  All  three  vary  in  time. 
Through  geometric  relations,  each  pixel's  location  is  deter¬ 
mined  by  its  time  of  observation. 

To  understand  these  relations,  a  general  examination  of 
the  earth,  satellite,  and  sensor  characteristics  is  first 
presented.  Second,  the  format  and  nature  of  the  pixel  data 
is  examined.  Third  the  relation  between  each  pixel  element 
and  its  corresponding  earth  location  is  derived.  Finally, 
an  algorithm  is  presented  to  convert  pixel  coordinates  to 
geographic  coordinates  and  vice  versa. 

Earth  Characteristics 

To  a  first  approximation,  the  earth  is  spherical  with 
the  same  volume  as  a  sphere  with  a  radius  of  6371  km 
[4:797].  Any  location  on  the  surface  of  the  earth  can  be 
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referenced  by  geocentric  latitude  and  longitude  coordinates 
with  respect  to  the  earth's  center  and  equatorial  plane. 

The  earth  rotates  from  west  to  east  about  its  north- 
south  axis  as  it  orbits  the  sun.  During  one  tropical  year 
(355.2422  days),  the  earth  makes  one  extra  revolution  of  360 
degrees  relative  to  the  stars  (366.2422  revolutions  in  all). 
The  mean  rate  of  revolution  is  360.9856  deg/day  (15.0411 
deg/hr  or  0.2507  deg/min)  [20:B6].  Viewed  from  a  fixed 
location  and  direction  in  space,  each  earth  location  will 
move  from  west  to  east  at  15.0411  deg  longitude/day  along 
its  parallel  of  latitude  (nearly  constant) . 

Actually,  the  earth  bulges  at  the  equator  and  sags  at 
the  poles  due  to  the  rotational  forces.  Its  shape  is  better 
modeled  as  an  oblate  spheroid  whose  radius  is  a  function  of 
latitude.  In  the  International  Astronomical  Union  model 
(lAU;  1976) ,  the  mean  equatorial  and  polar  radii  are 
6378.140  km  and  6356.755  km,  respectively  [20:K13].  This 
represents  an  eccentricity  of  0.08182.  At  NOAA,  radii 
values  of  6378.144  km  and  6356.759  km  are  used  [10:C-10]. 

Even  the  ellipsoid  model  is  an  approximation.  The 
equator  is  slightly  elliptical  with  an  eccentricity  of  about 
1  X  10'^  [12:25].  This,  combined  with  irregularities  in  the 
earth's  crust,  cause  the  local  acceleration  of  gravity  to 
vary  [12:38].  The  net  effect  is  to  produce  small  changes  in 
the  satellite  motion  and  viewing  geometry. 


Satellite  Orbits 


The  TIROS-N  series  of  satellites  are  in  polar  orbits, 
crossing  the  equator  twice  and  passing  over  both  polar 
regions  once  on  each  orbit.  The  rotation  of  the  earth  under 
the  satellite  causes  the  suborbital  track  to  shift  westward. 

The  gravitational  attraction  of  the  earth's  equatorial 
bulge  causes  the  line  of  nodes  of  the  orbital  plane  to 
process  eastward  (rotate  about  the  earth's  axis).  For 
sun-synchronous  satellites,  the  orbital  inclinations  are 
chosen  to  make  the  rate  of  precession  equal  to  the  change  in 
orientation  with  respect  to  the  sun  as  the  earth  proceeds  in 
its  orbit  about  the  sun.  On  average,  this  rate  is  0.9856474 
deg/day  (east  longitude)  [20:B6].  The  TIROS-N  satellites 
are  sun-synchronous,  passing  over  a  given  earth  location  a- 
about  the  same  local  time  every  day. 

The  designed  TIROS-N  constellation  consists  of  two 
satellites  in  sun-synchronous  orbits.  The  nominal  orbital 
parameters  of  height,  inclination,  period,  and  increment 
(longitude  shift  between  successive  orbits)  are  summarized 
in  Table  1. 

The  two  circular  orbits  intersect  at  right  angles,  so 
that  one  satellite  follows  the  other  about  six  hours  later. 
The  orbits  are  chosen  to  have  a  morning  satellite  that 
passes  overhead  on  a  descending  pass  between  0600  and  1000 
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Table  1.  TIROS-N  Nominal  Orbital  Parameters  [1:3-5] 


Heiaht 

(km) 

Inclination 

(deg) 

Period 

(min) 

Increment 

(deg  W) 

Orbits 

(#/day) 

833 

98.739 

101.38 

25.40 

14.18 

870 

98.899 

102.37 

25.59 

14.07 

local  mean  solar  time  (LMT) ,  and  an  afternoon  satellite  that 
passes  overhead  on  an  ascending  pass  between  1400  and  1800 
LMT  (see  Figure  1)  [1:3-6]. 

New  satellites  are  put  into  orbit  to  replace  aging  ones 
in  order  to  maintain  the  constellation.  The  replaced  satel¬ 
lites  may  continue  to  transmit  data  as  long  as  operational. 

As  of  this  writing,  NOAA-10  and  NOAA-11  are  the  primary 
satellites  of  the  TIROS-N  constellation.  NOAA-9  continues 
to  transmit  when  not  in  conflict  with  NOAA-11.  A  sample  of 
orbital  data  is  given  in  Table  2. 

Orbital  data  may  be  obtained  from  the  NOAA  Direct 
Readout  Users  Electronic  Bulletin  Board  (EBB)  [17:7-3]  or 
the  TBUS-1  and  TBUS-2  bulletins  transmitted  through  the 
National  Weather  Service  Communications  Center  (KWBC) 

[1:5-2]. 

Small  departures  from  circular  sun-synchronous  orbits 
occur  due  to  initial  orbit  injection  errors  and  subseguent 
perturbations.  Initial  orbital  parameters  are  expected  to 
meet  the  target  values  within  the  following  limits  [15:21]: 
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Figure  1.  TIROS-N  Constellation  Orbits  [1:3-7] 
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Table  2.  Sample  Orbital  Prediction  Data  [16] 


Orbital  Prediction  Data  for  August  1989 


Orbit  Number 

Equator  Crossing  Time  (Z) 
Long.  Ascending  Node  (deg) 
Nodal  Period  (minutes) 
Frequency  (MHz) 

Increment  (deg) 


NOAA-9 

NOAA-10 

NOAA-11 

23879 

14906 

4380 

0123.45 

0016.35 

0023.89 

129. 61W 

71.03W 

158. 76W 

102.0426 

101.2340 

102 . 1064 

137.62 

137 . 50 

137 . 62 

25.51 

25.31 

25.53 

Altitude  (average)  :  +  15  km 

Inclination  :  ±  0.150  deg 

Apogee/Perigee  Difference  ;  less  than  56  km 

Dominant  perturbations  include  the  gravitational  attraction 
of  the  sun  and  moon,  drag,  and  the  varying  gravity  of  an 
irregular  earth. 

The  spacecraft's  attitude  is  controlled  by  the  Attitude 
Determination  and  Control  System  (ADACS) .  The  ADACS  uses 
three  orthogonal  Reaction  Wheel  Assemblies  (RWA)  in  combina¬ 
tion  with  earth  and  sun  sensors  to  maintain  a  fixed  orienta¬ 
tion  with  respect  to  local  nadir  (straight  down  from  the 
satellite  towards  the  center  of  the  earth) .  This  orienta¬ 
tion  is  maintained  within  0.2  degrees. 

Departure  from  ideal  circular  orbits,  combined  with 
local  perturbations,  produce  small  attitude  errors.  Analys- 
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is  has  shown  that  attitude  is  generally  maintained  within 
0.12  degrees,  with  only  brief  perturbations  of  up  to  0.2 
degrees  [15:5].  Assuming  a  spherical  earth  with  radius  of 
6371  km,  these  correspond  to  horizontal  deviations  at  the 
subpoint  of  1.74-2.91  km  and  1.82-3.04  km  at  the  nominal 
satellite  altitudes  of  833  km  and  870  km,  respectively. 

Sensor  Characteristics 

The  current  TIROS-N  spacecraft,  called  Advanced  TIROS-N 
(ATN) ,  carry  several  sensor  systems  including  the  following: 
TIROS  Operational  Vertical  Sounder  (TOVS) ,  consisting  of  the 
High  Resolution  Infrared  Radiation  Sounder  (HIRS/2) ,  Strato¬ 
spheric  Sounding  Unit  (SSU) ,  and  Microwave  Sounding  Unit 
(MSU) ;  Advanced  Very  High  Resolution  Radiometer  (AVHRR) ; 
Space  Environmental  Monitor  (SEM) ;  Data  Collection  System 
(DCS);  and  Search  and  Rescue  System  (SAR)  [13:2-3].  Table  3 
summarizes  their  functions. 

The  earth  images  are  derived  from  the  AVHRR  sensor. 

This  is  the  primary  data  source  for  both  High-Resolution 
Picture  Transmission  (HRPT)  at  full  resolution,  and  Automat¬ 
ic  Picture  Transmission  (APT)  at  reduced  resolution. 

The  AVHRR  is  a  scanning  radiometer,  measuring  radiant 
energy  in  four  (AVHRR/1)  or  five  (AVHRR/2)  spectral  channels 
(see  Table  4)  [15:44].  NOAA-10  carries  the  last  four- 
channel  AVHRR,  while  NOAA-9  and  NOAA-11  have  five-channel 
AVHRR  [8:1-34,1-40,6-3;  10: Appendix  B  (NOAA-F/9  and 
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Table  3.  Advanced  TIROS-N  Primary  Sensors  [17:3-4] 


1.  TIROS  Operational  Vertical  Sounder  (TOVS) 

a.  Measures  temperature  profile  of  earth's 
atmosphere  from  surface  to  10  millibars 

b.  Measures  water  vapor,  ozone,  carbon  dioxide,  and 
oxygen  content  of  earth's  atmosphere 

2.  Advanced  Very  High  Resolution  Radiometer  (AVHRR) 

a.  Measures  radiation  in  visible  and  infrared 
spectra 

b.  Primary  source  of  data  for  High  Resolution 
Picture  Transmission  (HRPT)  and  Automatic  Picture 
Transmission  (APT)  data 

3.  Space  Environment  Monitor  (SEM) 

Detects  radiation  from  space  at  various  wavelengths 

4.  Data  Collection  System  (DCS) 

Collects  data  from  earth-based  environmental  plat¬ 
forms 

5.  Solar  Backscatter  Ultraviolet  Radiometer  (SBUV) 

a.  Measures  total  ozone  content  and  vertical 
distribution 

b.  Carried  on  spacecraft  in  afternoon  ascending 
orbits 

6.  Search  and  Rescue  (SAR) 

Detects  and  locates  emergency  beacon  radio  signals 
carried  on  ships  and  aircraft 

7.  Earth  Radiation  Budget  Experiment  (ERBE) 

Experiment  to  gather  data  on  earth's  average  radia¬ 
tion  budget 
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Table  4.  AVHRR  Spectral  Channels  [15:44] 


AVHRR/1  AVHRR/ 2 

(4  channels)  (5  channels) 


Channel 

Wavelenath 

Channel 

Wavelenath 

(microns) 

(microns) 

1 

0.55  - 

0.68 

1 

0.58  -  0.68 

2 

0.725  - 

1.10 

2 

0.725  -  1.10 

3 

3.55  - 

3.93 

3 

3.55  -  3.93 

4 

10.5  - 

11.5 

4 

10.3  -  11.3 

5 

(Channel  4) 

5 

11.5  -  12.5 

NOAA-G/10) 

,  2 -Rev,  3-] 

Rev]  (Note: 

NOAA-C/7 

and  NOAA-G/10  should 

be  listed 

as  having 

five-  and 

four-channel  AVHRR,  respec- 

tively,  as  in  Appendix  B) . 

The  Instantaneous  Field  of  View  (IFOV)  for  all  channels 
is  1.3  mr  +  0.1  (Imr  =  10'^  radians).  The  IFOV  of  all  chan¬ 
nels  are  also  coincident  within  O.l  mr  [15:44].  This  par¬ 
ticular  IFOV  was  chosen  so  that  subsequent  scans  are  con¬ 
tiguous  at  the  subpoint  as  the  satellite  orbits  the  earth 
[15:43] . 

The  AVHRR  utilizes  a  rotating  mirror  to  scan  through 
360  degrees  at  a  rate  of  360  revolutions  per  minute  (rpm) . 
Each  scan  is  perpendicular  to  the  satellite's  direction  of 
motion,  crossing  through  the  satellite  subpoint  (nadir) 
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[15:43].  The  scan  angle  represents  the  mirror  angle  measur¬ 
ed  in  reference  to  the  nadir  direction. 

Thus,  the  AVHRR  is  continuously  sensing  radiation  in 
all  channels  as  the  IFOV  rotates  through  360  degrees.  Since 
the  mirror  rotates  at  360  rpm,  this  produces  6  lines  per 
second.  The  scan  direction  is  from  the  space-view  side  of 
the  satellite  across  the  earth  toward  the  sun  side  (from 
right  to  left)  [15:46]. 

One  scan  across  the  earth  produces  a  line  of  data 
referred  to  as  a  "scan  line."  Each  scan  line  is  skewed  by 
the  forward  motion  of  the  satellite  during  the  scan. 

The  sensor  actually  scans  along  a  scan  vector  which  is 
perpendicular  to  both  the  satellite's  position  and  velocity 
vectors.  Projected  onto  the  earth's  surface,  this  forms  a 
great  circle  arc  at  right  angles  to  the  satellite  track  (see 
Figure  2) . 

Applying  the  four-parts  formula  (cosacosC  =  sinacotb  - 
sinCcotB)  to  Figure  2,  the  scan  skew  7  is  given  by 

tan7  =  tan^/sinCa^ )  (1) 

where 

♦  =  earth  center  angle  (radians) 

=  satellite  angular  progress  (radians) 

All  distances  along  the  scan  line  are  proportional  to  the 
ECA  by  a  factor  of  sin7.  This  is  used  in  calculating  the 
along-line  resolution. 
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Legend 

Points  Angle 

E:  Earth  location  7;  Scan  skew 

S':  Satellite  subpoint 

Arcs 

Earth  Center  Angle  (ECA)  from  S '  to  E 
he:  Angular  distance  along  satellite  orbit 


Figure  2.  AVHRR  Scan  Line  Skew 


During  half  a  scan  line  after  nadir,  the  satellite 
travels  an  angular  distance  he  =(27r/P)t,^2  '  where  P  is  the 
satellite  period,  and  t^^j  is  the  time  to  complete  half  a 
scan  line  (0.02563  seconds).  The  ECA  is  also  at  a  maximum 
'^max  •  nominal  870  km  orbit  with  a  period  of  102.37 


min,  =  13.86  deg  (0.2420  rad)  and  ^6  =  0.001502  deg 

(1.5886  mr) ,  giving  a  scan  skew  of  89.63  deg. 

Some  authors  [3;  6]  have  referenced  the  scan  angle  to 
the  satellite  subpoint  at  the  center  of  the  scan  line. 
Actually,  the  subpoint  continues  to  move  during  the  scan. 

The  skewed  nature  of  the  scan  line  produces  an  apparent  scan 
angle  that  has  more  to  do  with  the  problem  of  displaying  the 
data  than  the  actual  operation  of  the  scan  mirror.  In  this 
work  the  scan  angle  is  defined  as  the  actual  mirror  angle 
measured  in  reference  to  the  subpoint  direction  at  the  time 
of  viewing. 

Sensor  Data 

The  TIROS-N  data  handling  subsystem  has  two  information 
processors:  the  TIROS  Information  Processor  (TIP)  for  low 

data  rates,  and  the  Manipulated  Information  Rate  Processor 
(MIRP)  for  high  data  rates.  The  data  flow  is  outlined  in 
Figure  3. 

The  TIP  formats  (and  controls)  the  data  output  from  the 
low-bit-rate  instruments  (HIRS/2,  MSU,  etc.)  and  the  teleme¬ 
try  data.  It  also  adds  a  time  code  accurate  to  one  mil¬ 
lisecond  (see  Figure  4) ,  identification,  and  synchronization 
[15:78] . 

The  MIRP  processes  and  formats  (and  controls)  data  from 
the  AVHRR,  adding  synchronization,  telemetry,  and  a  time 
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Figure  3.  TIROS-N  Data  Flow  Diagram  [15:7] 


Figure  4.  TIP  Time  Code  Format  [15:9] 


code.  TIP  data  is  also  included  (except  for  APT  process¬ 
ing)  .  MIRP  output  is  used  for  the  following:  real-time 
HRPT  and  APT  transmissions,  tape-recorded  Global  Area  Cover- 
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age  (GAC) ,  and  Local  Area  Coverage  (LAC)  [15:9-11].  Table  5 
summarizes  the  MIRP  output  characteristics. 


Table 

5 .  MIRP  Output 

Characteristics 

[15:11] 

HRPT/LAC 

GAC 

APT 

Form  of 
data 

Serial  digital  bit  stream  of 
10-bit  words,  MSB  first 

Analog  AM  on 
2,400  Hz 
subcarrier 

Line  rate 

6/sec 

2/sec 

2/sec 

Word  rate 
(words/sec) 

66,540 

6,654 

4,160  before 
digital-to- 
analog 
conversion 

Number  of 

AVHRR 

channels 

s 

5 

2 

Words  of 
earth  data 
per  line 
per  channel 

2048 

409 

909  before 
digital-to- 
analog 
conversion 

Processing 
of  AVHRR 
data 

Formatting 

only 

Resolution 

reduction; 

formatting 

Resolution 

reduction; 

formatting 

Other  data 
than  AVHRR 

TIP;  time  code 

TIP;  time  code 

Minute  marks 

Calibration 

wedge 

When  the  AVHRR  scan  mirror  reaches  a  predetermined 
pre-earth  position,  a  synchronization  pulse  is  generated  and 
sent  to  the  MIRP.  This  event  becomes  the  reference  for  all 
subsequent  event  timing  [15:12]. 


Event  timing  is  based  on  a  0.9984  megahertz  (MHz) 
internal  clock  [14].  At  this  frequency,  the  clock  makes 
0.9984  X  10 counts  every  second,  or  one  count  every 
1.001602564  X  10*  seconds.  Events  may  also  be  referenced 
in  Standard  Time  Units  (STU) ,  where  1  STU  =  1000  counts. 
Table  6  summarizes  the  timing  of  major  events. 


Table  6.  MIRP  Event  Timing  [14] 


Count 

STU 

Event 

0 

0 

Frame  sync,  bit  1,  word  1 

8,600 

8.6 

First  pulse  of  earth  data 

34,200 

34.2 

Nadir 

59 , 775 

59 .775 

Last  pulse  of  earth  data 

68,400 

68 . 4 

End 

The  MIRP  issues  sampling  commands  to  the  AVHRR,  making 
instantaneous  measurements  of  the  instrument  readings.  Each 
reading  is  coded  as  a  ten-bit  word  with  the  most  significant 
bit  (MSB)  first.  In  each  scan,  the  AVHRR  produces  2,048 
samples  of  earth  data  for  each  channel,  along  with  other 
reference  samples  as  indicated  in  Figure  5  [15:10-13]. 

As  illustrated  in  Figure  6,  in  each  AVHRR  scan  line, 
2047  words  (2048  -  1)  are  produced  in  51,175  time  clock 
counts  (59,775  -  8,600),  or  one  word  every  25  counts 
(25.04006  X  10’*  seconds).  The  mirror  rotates  at  6  rev/sec. 
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Figure  5.  AVHRR  Data  Sampling  Areas  [15:13] 


-IFOV- 


IFOV- 


— word- 


I — word - j 

Word 

2  3 


8600  8625  8650  8675  8700 

Count 


2047 

H 

59775 


Figure  6.  AVHRR  Scan  Line  IFOV 
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SO  there  are  0.9439882  x  10  radians  (mr)  per  word,  or 
110.7692308  degrees  per  line  of  2048  words.  Since  the  IFOV 
is  l.j  mr,  the  views  of  consecutive  pixels  (words)  overlap 
by  about  0.36  mr  (27%  of  the  IFOV). 

HRPT/^PT  Processing 

The  MIRP  also  performs  the  function  of  processing  two 
channels  of  AVHRR  data  into  a  reduced  resolution  format  (and 
lower  data  rate)  for  APT  transmission.  Only  one  AVHRR  line 
in  three  is  processed  for  each  of  the  two  APT  channels, 
producing  two  lines/sec  of  APT  data.  The  APT  samples  are 
then  converted  to  an  analog  format  with  an  accuracy  of  8 
MSBs  of  the  original  10-bit  AVHRR  words  [15:63].  The  phas¬ 
ing  of  the  MIRP  output  is  summarized  in  Figure  7, 

The  data  reduction  algorithm  divides  each  scan  into 
five  regions  according  to  the  angular  distance  from  nadir  as 
shown  in  Table  7.  This  algorithm  reduces  2,048  HRPT  words 
to  909  APT  words  per  channel,  providing  a  nearly  uniform 
resolution  of  4  km  along  the  scan  line. 

The  correspondence  between  APT  and  HRPT  pixels  is  given 
in  Table  8.  The  scan  angle  is  referenced  to  nadir  (0  deg). 
Pre-nadir  angles  are  defined  as  negative  and  post-nadir 
angles  defined  as  positive. 
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Figure  7.  MIRP  Output  Phasing  [15:14] 


Table  7.  HRPT/APT  Data  Reduction  Algorithm  [15:16] 
Region  Scan  Angle  *  Reduction  Method 


degrees  (+) 

1 

0.0  -  16.9 

Average  4  samples 
(contiguous) 

2 

16.9  -  34.8 

Average  2  samples, 
skip  1 

3 

34.8  -  43.8 

Average  2  samples 
(contiguous) 

4 

43.8  -  48.8 

Average  1  1/2  samples 
(A+B)/2;  (B+C)/2 

5 

48.8  -  55.4 

Retain  original 
(No  reduction) 

Measured 

from  nadir 
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Table  8.  HRPT/APT  Pixel  Correspondence 


HRPT  APT 


Scan  Ancle 

Words 

Pixel 

Ratio 

Pixel 

Words 

-55.38 

to 

-48.84 

121 

0-120 

1:1 

0-120 

121 

-48.84 

to 

-43.81 

93 

121-213 

3:2 

121-182 

62 

-43.81 

to 

-34.83 

166 

214-379 

2:1 

183-265 

83 

-34.83 

to 

-16.98 

330 

380-709 

3:1 

266-375 

110 

-16.98 

to 

16.98 

628 

710-1337 

4:1 

376-532 

157 

16.98 

to 

34.83 

330 

1338-1667 

3:1 

533-642 

110 

34.83 

to 

43.81 

166 

1668-1833 

2:1 

643-725 

83 

43.81 

to 

48.84 

93 

1834-1926 

3:2 

726-787 

62 

48.84 

to 

55.38 

121 

1927-2047 

1:1 

788-908 

121 

Total 

2048 

0-2047 

(4km) 

0-908 

909 

Data  Format 

HRPT  is  a  real-time,  high-data-rate  digital  transmis¬ 
sion  (see  Table  9  and  Table  10) .  The  HRPT  transmission  con¬ 
tains  all  spacecraft  instrument  data  [10: 10-rev].  This  con¬ 
sists  of  output  from  the  five  AVHRR  channels,  multiplexed 
with  the  TIP  output  (including  telemetry). 

The  HRPT  line  and  minor  frame  formats  are  summarized  in 
Figure  8  and  Table  11.  Three  minor  frames  make  up  a  major 
frame.  The  TIP  data  is  updated  only  once  every  major  frame 
(every  three  minor  frames)  [10:18,20-23;  15:71-75]. 


Table  9.  HRPT  Characteristics  [10:13] 

Line  rate  6  lines/second 

Carrier  modulation  Digital  split  phase, 

phase  modulated 

Transmit  frequency  1698.0  MHz  or  1707.0  MHz 

Transmit  power  8  watts  nominal 

EIRP  (approximate)  39.0  dBm 

Polarization  Right  hand  circular 

Spectrum  bandwidth  2.4  MHz  (3dB) 

Table  10.  HRPT  Parameters  [10:17] 


Manor  Frame 

Rate 

2  frames/sec 

Format 

3  minor  frames 

Minor  frame 

Rate 

6  frames/sec 

Format 

See  Figure  8 

Word? 

11,090 

Word 

Rate 

66,540  words/sec 

Words/ channel 

2 , 048/line 

Bits 

10 

Order 

Bit  1  =  MSB 

(Most  Significant  Bit) 

Bit  10  =  LSB 

(Least  Significant  Bit) 

Bit 

Rate 

665,400  bits/sec 

Format 

Split  phase 

VAC<  PAtA  MMMM 

AACIKAMOAfA  AVHAAONLV 

UltMttPV  OAtA  ItMDIlM 


-  iO(«riricATtoii  a«OAM 


Hi  MM«OII  fAAMf  LfMOtH  t«.M»«OAOI 
Oi  iHAIIMlMOAfAAMCIHAMAJOAfAAMt 
Qt  MMOA  f  AAM(  AAlf  •  rAAMMCOMO 
M)  WOAOLfMGIH  MilltAMAO 
m  ALL  VAAilAAC  MTM  OCOMt  CMt  •AAt 


TIMNOAO  AllOCAItOM 


l«  AAMT  CAlMAAf  M>N 

•  CHANMCl S lAAOf T 
icuroff  iUACOMI 

I  CMAMMIlAlAAGiT 
IIMT  (A  AT  CUaOM) 

•  CHAMNCL I TAAOil 
lik*  (A  AT  AUlCOU) 

•  aUNNfLIAATCll 
TIMA 

M  aaari 


1  aWMIDIIIAlLOCAIIOM 

MtioaMo 

tmmmono 

1  tvaicio 

IS  MAMA  ID 

4  >  tAACtCAArt  AMNIIU 

•  MtVMCMAilKIII 

•  AATAA 

M  OAtA 1 

KAAAil 

Figure  8.  HRPT  Line  Format  [10:18] 


APT  is  a  realtime,  mediuro-data-rate  analog  transmis¬ 
sion  ^see  Table  12  and  Table  13).  The  APT  transmission  con- 
contains  MIRP  processed  data  from  two  ground  selectable 
AVHRR  channels,  labeled  channels  A  and  B  [10: 10-rev]. 

Channels  A  and  B  are  multiplexed  with  telemetry  and 
space  data,  then  converted  to  an  analog  signal  which  ampli¬ 
tude  modulates  (AM)  a  2400  Hz  subcarrier  for  analog  trans¬ 
mission  on  a  frequency  modulated  (FM)  carrier  [10: 10-Rev; 
15:16]  . 
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Table 


HRPT  Minor  Frame  Format  [10:20] 


Funclloo  I  No.  of  Words 


Word 

Posltlos  1 


Bit  No. 
4  S  6  T 


101000010  0 
0  10  110  111  1 
110  10  1110  0 
0  110  0  1110  1 
lOOOOOlll  1 
001001010  1 


Plus  word  coda  !■  meaning 


First  60  Uta  from  a  63-bltPNU) 
genaralor  started  In  the  all  I'a 
state.  The  generator  poly- 
nomlnal  la  f  xS  X  +  1 


Bit  It  0  -  Internal  synoi  1  ■  AVIllUl  ayno 
Bite  3  li  S|  00  ■  not  usedf  01  -  minor  frame  1| 

10  >  minor  frame  >,  11  •  minor  frame  3 
Bite  4-7t  spacecraft  addreasi  bit  4  •  MSB,  bit  7  °  L8B 
Bit  St  0  •  frame  atablai  1  >  frame  resync  occurred 
DIU  0-10|  aparei  bit  0-0,  bit  10-1 
Spare  wordi  bit  symbols  undefined 


Bite  1-Ot  binary  day  counit  bit  1  -  MSBt  bit  0  -  1.SB 

Bit  10|  0|  spare 

Bite  l-3t  all  O'at  spare  1,.  0,  1 

Bits  4-lOt  pert  of  binary  mauc  of  day  eounli  bit  4  -  MSB 
of  meeb  count 

Bit  1-lOt  part  of  binary  msec  of  day  ccunlt 
Bit  1-lOt  remainder  of  binary  msec  of  day  eaunl; 


Ramp  calibration  AVilRR  channel  1 
Ramp  calibration  AVIIRR  channel  3 
Ramp  calibration  AVIIRR  channel  3 
Ramp  callbrallcn  AVIIRR  channel  4 


tomperalure  ^ 

V»  rAianaiM  vaIua 
AVHRB  patcil  temperaluro  reieranoe  value 

000  0  00000  1  apare 


10  worda  of  latenwl  target  data  from  each  AVHRR  cb  3, 
4,  and  6.  Tbeae  data  are  time  multipleaed  aa  ch  3 
(word  1).  eh  4  (word  1),  oh  6  (word  1),  oh  3 
(word  31,  ch  4  (word  3),  ch  6  (word  1),  ate. 


10  worda  of  spaee-aoaa  data  from  each  AVHRR  channel 
1,  3,  3,  4,  ami  0.  Theae  data  are  lima  multipleaed  aa 
oh  1  (word  1),  eh  3  famrd  1),  oh  3  (word  1),  ch  4 
(word  1),  oh  3  (word  1),  ch  1  (word  3),  oh  3  (word  3), 
ch-  3  (word  3),  ch  4  (word  3),  oh  6  (word  3),  etc. 


Bit  l(  0  -  AVHRR  sync  eaiiyi  1  -  AVHRR  ayno  lata 
Bila  3-10|  8-blt  binary  count  of  0.  »B84-Mlls  periodai 
bit  3  -  MSB,  bit  10  -  LSB 


(2|  Aa'meaaurod  by  a  platinum  reaislanoe  thermometer  embedded  In  the  housing. 


Table  11  (continued) 


Function 

No.  of  Worda 

Word 

Poaillon 

UU  No. 

ia34S«7a9I0 

Plua  Word  Coda  k  Maaninf 

Tip  (JaU 

S20 

104 

llie  630  ivonia  eonUtn  flv«  frames  of  TIP  data  (104  TIP  I 

data  worda/frama) 

BIta  I'Oi  aaact  format  aa  (eoeralad  by  TIP 

Bit  9i  avea  parlQi  aback  ovar  bita  1-8 

623 

Bit  lOi  -  bit  1 

Spare  wordn 

624 

101  00  0111  ON 

127 

111  0  0  0  1  0  1  1 

Darivad  by  invarllnc  tha  output 

000010111  1 

of  a  102S-blt  PN  aequenca  pro- 

tl9*9 

1AIIAAA1I  1 

vidad  by  a  faadbaok  ahift  ragla- 

ilOlOlOOl  0 

1 

tar  ganaratlM  tha  polynomlnali 
’  X**  ♦  X“  X*  +  X  +  1 

1 

748 

1 

100-101101  0 

Tha  caaaralor  la  klartad  In  tha 

749 

iiooiobol  0 

I'a  alata  at  tha  beginning  of 

760 

100000000  0^ 

word  7  of  each  minor  frame. 

Earth  data 
(AVIIRR) 

10.240 

I 

Ck  1  -  Bampla  1  '< 

Ch  2  -  Rampla  1 

Ch  3  -  Sam^a  1 

Ch  4  -  Bampla  1 

Ch  S-bmplal 

Ch  1  -  Sample  2 

* 

Ch  S- Bampla  2047 

Ch  1  -  Sample  2048 

Ch  2 -Bampla  2048 

Ch  3 -Sample  2048 

Ch  4  -  temple  2048 

Ch  8  -  temple  2048  > 

Each  minor  frame  eontalna  the 
data  obtained  during  one  earth 
acaa  of  the  AVHRIt  aanaor. 

^  Tha  data  from  tha  five  aanaor 
chauwla  of  tha  AVIIRR  are 
lima  iiiulllplaxad  aa  Indicated 

Auxiliary 

sync 

100 

10.991 

10.992 

10.993 
10,094 

1 

11,089 

11,090 

1  1.  1  1  1  0  0  0  1  o' 

I  1  1  1  1  1  0  0  1  1 

0  110-11010  1 

1  0  1  0  1  1  1  1  0  1 

0  1  1  t  1  1  0  0  0  0 

1  1  1  1  0  0  1  1  0  0^ 

Derived  from  tha  .nonlnverted 
output  of  a  1023-bit  PN  aa- 
quenca  provided  by  a  feedback 
ahin  regialar  ganermting  the 
polynomlnal:  X  X^  «  X^  •  X  •  1 

Tha  generator  la  alerted  in  the 
all  I'a  alata  at  the  beginning  of 
word  10,891 

As  of  this  writing,  the  primary  spacecraft  (NOAA-10  and 
NOAA-11)  are  transmitting  one  visible  and  one  infrared  (IR) 
channel  during  spacecraft  day,  and  two  IR  channels  during 
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Table  12,  APT  Characteristics  [10:11] 


Line  rate 

2  lines/sec 

Data  resolution 

4  km  (nearly  uniform) 

Carrier  modulation 

Analog 

Transmit  frequency 

137.50  MHz  or  137.62  MHz 

Transmit  power 

6  watts  nominal 

Polarization 

Right  hand  circular 

Subcarrier  frequency 

2.4  kHz 

Carrier  deviation 

±  17  kHz 

Low  pass  filter 
(ground  station) 

1.4  kHz  7th-order  linear 
recommended 

Synchronization 

Channel  A:  7  pulses  at  1040 
pps  with  50%  duty  cycle 

Channel  B:  7  pulses  at  832 
pps  with  60%  duty  cycle 

spacecraft  night  via  APT.  NOAA-9  transmits  one  visible 
channel  and  one  IR  channel  throughout  its  orbit  [16], 

The  APT  line  and  frame  formats  are  summarized  in 
Figure  9,  Figure  10,  and  Figure  11.  Note  that  one  frame 
consists  of  128  lines.  Each  APT  frame  contains  one  teleme¬ 
try  frame  of  16  wedges  repeated  on  8  consecutive  lines,  both 
after  channel  A  and  after  channel  B  data. 

Minute  markers  are  also  included  in  four  consecutive 
lines  (two  white,  two  black) .  While  the  time  of  each  mark 
is  not  specified,  limited  observation  found  the  start  of 
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Table  13.  APT 

Parameters  [10:13] 

Frame 

Rate 

1  frame  per  64  seconds 

Format 

128  lines 

Line 

Rate 

2  lines/sec 

Format 

See  Figures  8,  9  and  10 

Words 

2,080 

Channels 

2  (selected) 

Word 

Rate 

4,160  words/sec 

Words/channel 

909 

Analog-to- 

Digital 

Conversion 

Accuracy 

The  8  Most  Significant  Bits 
(MSBs)  of  the  10-bit  AVHRR 
words 

Low-oass 

filter 

Type 

3rd-order  Butterworth- 
Thompson 

Bandwidth 

2400  Hz  (3  dB) 

each  minute  indistinguishable  from  the  first  white  line  in 
the  four-line  sequence.  This  work  assumes  that  the  beginn¬ 
ing  of  each  minute  coincides  with  the  start  of  the 
synchronization  pulse  of  channel  A  in  each  line  with  the 
first  white  minute  marker. 

To  process  the  APT  data  on  a  digital  computer,  the 
analog  APT  signal  received  must  be  reconverted  to  digital 
form  after  demodulation.  This  is  accomplished  by  sampling 
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Stiected  ChiniMl  A 
Video  Oeti 

900  Words 


Spece  Oeti/Mlnute  Herktrs-A7  Words 


Sync  A- 39  Words 

TtItMtry  DeU-4S  Words. 


1 


Selected  Chennol  B 
Video  Oeu 

909  Words 


Speco  OeU/NinuU  Nirkers-47  Words 

Syec  1-39  Words 

Teleaotry  DeU-4S  Words _ 


1.  Equivilent  ootput  dloltal  deU  rate  Is  4160  Hords/sacond 

2.  Video  line  rata  -  2  llnet/socood 

3.  APT  fraiH  slia  -  128  lines 

4.  Any  two  Of  the  five  AVHM  channels  nay  be  selected  for  use 
$.  Sync  A  Is  a  1040-Hi  square  wave  -  7  cycles 

6.  Syne  8  is  a  aia-pps  pulse  train  -  7  pulses 

7.  each  of  16  telcMtry  points  are  repeated  on  B  successive  lines 

8.  Minute  aurkers  are  repeated  on  4  succassiva  lines,  with  2  lines 
black  and  2  lines  white 


Figure  9.  APT  Line  Format  [10:14] 


Af T lYNC  A 
41  cycilt  W  IMi  Hit 


AfTSVMCI 

<7PULttl«i3lppM 


unjinjuiRTi _ i 

I  1 

■Lririririjnja^ 


MAXMI*f?«n 


HAXMl-tVt  W 


•  4  t  tl  1« 


M  »  »  M 


m  rrne  a  muemt  chanml-a  data 
(»  fVNC  •  miCfOU  CHAMMl-t  DATA 


Figure  10.  APT  Synchronization  Foirmat  [10:16] 
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Figure  11.  APT  Frame  Format  [10:15] 


the  analog  signal  at  regular  intervals  and  assigning  a  digi¬ 
tal  value  to  the  amplitude  measured.  Comparing  digital  HRPT 
data  and  digitized  APT  data  has  shown  good  correlation  [21]. 
It  should  be  noted  that  the  digital  information  derived  is 
an  approximate  representation  of  the  original  APT  digital 
data  before  conversion  to  analog.  The  best  possible  case 


would  be  synchronized  sampling  at  the  APT  data  rate 
(4160  Hz)  to  recover  the  8  MSBs  of  the  original  10-bit 
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words.  A  sample  rate  higher  (or  lower)  than  the  data  rate 
has  the  effect  of  interpolating  (or  smoothing)  the  data. 

For  this  study,  a  sampling  rate  of  9600  Hz  was  used. 

The  resulting  4800  samples  per  line  represent  an  interpola¬ 
tion  of  the  original  2080  words  per  line.  At  this  sampling 
rate,  the  ratio  of  samples  to  words  is  30:13  (2.3077:1). 

Earth/Satellite  Viewing  Geometry 

The  basic  earth/satellite  viewing  geometry  lies  in  a 
plane  determined  by  the  satellite  position  S,  earth  location 
E,  and  earth  center  0  (See  Figure  12).  The  off-nadir  view¬ 
ing  angle  6,  earth  center  angle  (ECA)  4',  local  zenith  angle 
Z,  satellite  radius  r,  and  range  distance  d  from  satellite  S 
to  earth  location  E  are  related  as  shown  in  the  following 
discussion . 

Applying  the  cosine  law  in  triangle  SOE  gives 

2  =  (;j2  +  ]-2  _  2drcos5  (2) 

where  is  the  earth  radius  at  location  E. 

Solving  for  d  and  using  the  negative  root  (which  is  the 
closer  distance)  relates  the  range,  d,  to  off-nadir  angle  <5: 

d  =  rcos<S  -  (Rj  ^  sin^5)’^^  (3) 

where  r  is  the  satellite  radius  (from  0)  and  Rg  is  the 
earth  radius  at  location  E. 

The  Zenith  angle,  Z,  is  an  exterior  angle  of  triangle 
SOE,  so  that 

Z  =  <5  +  'K  (4) 
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Legend; 


Points 

Distances 

S: 

Satellite 

h; 

Satellite  height 

S '  : 

Subpoint 

d: 

Range  from  S  to  E 

E; 

Earth  Location 

R; 

Earth  radius 

0: 

Earth  Center 

r: 

Satellite  radius 

H: 

Earth  Horizon 

(r  =  R  +  h) 

Angles 

6:  Off-nadir  angle 

Earth  Center  Angle  (EGA) 
Z:  Zenith  angle 


Figure  12.  Earth/Satellite  Viewing  Geometry 
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Applying  the  law  of  sines  in  plane  triangle  SOE  also  yields 

sinZ  =  (r/Rg  )siniS  (5) 

Substituting  Equation  4  into  Equation  5  gives 

sin((5  +  4')  =  (r/R£)sin<S  (6) 

Using  the  identity  sin(A+B)  =  sinAcosB  +  cosAsinB  in 
Equation  6  gives  5  in  terms  of 

tan<5  =  sin't/[(r/R£  )  -  cos'i')  ]  (7) 

Combining  Equations  3  and  6  relates  d  to  5  and  $: 

d  =  rcos<S  -  RjCOs((S  +  4')  (8) 

These  general  relations  are  valid  regardless  of  the 
shape  of  the  earth  or  the  path  of  the  satellite.  The  fol¬ 
lowing  sections  will  apply  these  relations. 

Reception.  Coverage,  and  Resolution 

The  relations  above  make  no  assumptions  about  the 
earth's  shape  or  the  satellite  orbit.  To  simplify  the 
discussion  of  satellite  reception,  coverage,  and  resolution, 
in  this  section  the  earth  is  assumed  to  be  a  sphere  with 
radius  R  =  6371  km  and  the  satellite  orbits  are  circular. 

These  assumptions  are  accurate  to  a  first  approxima¬ 
tion.  The  earth  radius  only  varies  from  about  6378  km  at 
the  equator  to  6357  km  at  the  poles  [10:C-10].  This  dif¬ 
ference  is  within  +15  km  (0.235%)  of  the  mean.  The  TIROS-N 
satellite  heights  vary  by  less  than  56  km  (6.5%)  [15:21]. 

The  satellite  track  and  footprint  are  illustrated  in 
figure  13.  The  subpoint  latitude  and  longitude  are  related 


41 


,  -  equator 


satellite 
track 


Legend 


Points 


Angles 


Earth  location  i: 

Satellite  subpoint  i': 

Ascending  node  at 
static  longitude 
North  pole  (rotational) 


Orbital  inclination 
Latitude  crossing 
angle  at  subpoint  S' 


Arcs 


♦;  latitude 

X ' :  static  longitude  from  N 
Earth  Center  Angle  (ECA) 

9:  Great  circle  arc  from  N  to  S' 

L:  Footprint  departure  (small  circle  arc) 


Figure  13.  Satellite  Track  and  Footprint 


to  the  angular  distance  traveled  by  the  satellite  after 
equator  crossing  as  follows: 

sin*5,  =  sins  sin(i)  (9) 

sinA'j,  =  cos^/cos*  (10) 

where 

=  subpoint  latitude 

A '5,  =  static  subpoint  longitude  measured  from  the 
ascending  node 

S  =  angular  distance  from  the  ascending  node  to  the 
satellite  subpoint 
i  =  orbit  inclination 

The  static  longitude  is  measured  from  the  fixed  posi¬ 
tion  of  the  ascending  node.  This  provides  a  convenient 
reference  for  geometric  purposes.  Since  the  earth  rotates, 
the  static  longitude  of  any  point  on  the  surface  changes 
with  time. 

The  satellite  crosses  the  parallel  of  latitude  at 
an  angle  i ' : 

cos(i')  =  sinA’/sin^  (11) 

The  departure  is  the  small  circle  distance  between  two 
points  along  a  given  parallel  of  latitude.  The  departure  L 
at  latitude  ♦  can  be  approximated  with  a  great  circle  arc  by 

sin(L/2)  =  sin^-^^ /sin(i' )  (12) 

where 

“  maximum  earth  center  angle 
1'  =  latitude  crossing  angle 

This  is  used  below  to  calculate  satellite  coverage. 
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Reception.  Since  APT  reception  is  essentially  line- 
of-sight,  the  satellite  is  only  visible  to  a  ground  station 
when  above  the  local  horizon  [17:7-6].  In  other  words, 
reception  begins  when  the  zenith  angle  decreases  to  90 
degrees  and  the  satellite  rises.  This  situation  is  il¬ 
lustrated  as  triangle  SOH  in  Figure  12 . 

Denoting  horizon  values  by  the  subscript  H,  the  earth 
center  angle,  ,  is  given  by  cos4'„  =  /r  .  Assuming  a 
spherical  earth  with  mean  radius  of  6371  km  and  nominal 
satellite  heights  of  833  km  and  870  km,  this  gives  EGA 
values  of  27.80  deg  (0.4857  rad)  and  28.38  deg  (0.4952  rad), 
respectively.  Using  Equation  3,  these  correspond  to  off- 
nadir  6^  values  of  62.18  deg  and  61.62  deg. 

The  great  circle  distance  from  earth  location  H  to 
subpoint  S*  is  R'P„  ,  where  is  in  radians.  Using  the 
values  previously  calculated  gives  horizon  distances  of 
3,094  km  and  3,155  km  for  the  nominal  satellite  heights  of 
833  km  and  870  km,  respectively. 

So,  the  satellite  is  visible  whenever  the  suborbital 
track  lies  within  a  radius  of  about  3,094-3,155  km  from  the 
earth  station.  In  general,  the  higher  the  satellite,  the 
greater  the  visibility. 

Coverage.  The  coverage  of  a  given  earth  location  (the 
frequency  with  which  it  is  visible  to  the  satellite  sensors) 
depends  on  the  swath  width  of  each  earth  scan  and  the  dis- 
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tance  the  track  shifts  on  each  orbit.  Swath  width  is  deter¬ 
mined  by  the  satellite  height  and  maximum  scan  angle,  while 
track  displacement  is  determined  by  the  orbital  inclination 
and  period. 

As  noted  previously,  the  maximum  AVHRR  scan  angle  is 
±  55.4  deg  from  nadir.  Using  Equation  4  with  the  nominal 
satellite  altitudes  of  833  km  and  870  km  gives  values  of 
±  13.15  deg  and  +  13.92  deg,  respectively.  These  correspond 
to  earth  surface  distances  of  ±  1,463  km  and  +  1,547  km  from 
the  satellite  subpoint  (swath  widths  of  2,926  km  and  3,094 
km)  . 

Recall  that  the  satellite  is  visible  at  an  earth  sta¬ 
tion  at  distances  up  to  3,094-3,155  km  from  subpoint.  This 
means  that  the  satellite  may  be  visible  to  an  earth  station 
without  providing  sensor  coverage  of  that  location. 

Since  data  can  be  received  when  the  satellite  is  within 
a  circle  surrounging  the  station,  maximum  coverage  is  avail¬ 
able  to  stations  on  the  satellite  track,  as  illustrated  in 
Figure  14.  This  provides  coverage  of  an  area  about  3,000  km 
X  6,250  km  (18,750,000  km^ )  centered  at  the  earth  station. 

The  longitude  distance  spanned  by  each  swath  at  the 
equator  is  slightly  greater  than  the  swath  width  due  to  the 
orbit  inclination  (see  Figure  13) .  Using  a  nominal  inclina¬ 
tion  of  98.8  deg  and  swath  widths  of  2,926  km  and  3,094  km, 
solving  Equation  13  gives  longitude  widths  at  the  equator 
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Figure  14 .  Maximum  AVHRR  Sensor  Coverage 


of  26.31  deg  and  28.09  deg,  or  2,960  km  and  3,134  km,  respe¬ 
ctively. 

As  discussed  earlier  (see  Earth  Characteristics) ,  the 
eastern  rotation  of  the  earth  causes  the  satellite  track  to 
shift  westward  at  the  rate  of  15.0411  deg/hr.  For  a  nominal 
satellite  period  of  102.37  min,  this  amounts  to  an  increment 
of  25.59  degrees  of  longitude  per  orbit.  With  a  swath 
spanning  about  28.09  deg  longitude  at  the  equator,  each 
earth  location  is  covered  at  least  twice  each  day  (one 
ascending  pass  and  one  descending  pass) . 

The  departure  (distance  between  degrees  of  longitude 
along  constant  latitude)  varies  with  latitude  ♦  in  propor- 
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tion  to  cost  [12:39].  Thus,  the  higher  the  latitude,  the 
less  the  westward  shift  due  to  rotation. 

The  longitude  difference  aa  spanned  by  a  given  depar¬ 
ture  L  is 

AA  =  L/COS*  (13) 

On  average,  the  number  of  passes  per  day  that  cover  an  earth 
location  at  latitude  ♦  is 

(passes/day)  =  (orbits/day) aa/180  (14) 

Solving  Equations  12,  13,  and  14  gives  the  satellite  cover¬ 
age  (passes/day)  for  both  satellite  reception  and  sensor 
coverage  using  the  appropriate  maximum  EGA  values  )  . 

The  average  number  of  reception  and  coverage  passes  per 
day  are  given  in  Table  14  at  various  latitudes  for  a  nominal 
870  km  orbit.  These  are  approximate  values  only.  The 
results  of  a  computer  simulation  over  a  30  day  period  pro¬ 
duced  similar  data,  though  deviations  up  to  1  pass/day  were 
noted  at  higher  latitudes. 

Resolution.  The  resolution  of  a  pixel  is  the  size  of 
an  area  on  the  earth's  surface  that  just  fills  the  IFOV  of 
the  sensor.  This  depends  on  the  angular  IFOV,  range  from 
satellite  to  the  area  imaged,  and  image  skew  due  to  off- 
nadir  viewino  »nd  the  curvature  of  the  earth. 

The  IFOV  can  be  expressed  as  the  difference  in  scan 
angle  between  the  beginning  and  the  end  of  a  pixel. 
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Table  14.  Earth  Station  Daily  Coverage  by  Latitude 


Latitude 

Coveraae 

Reception 

degrees  (N/S) 

passes/ day 

passes/day 

0 

2.2 

4.5 

15 

2.3 

4.7 

30 

2.5 

5.2 

40 

2.9 

5.9 

50 

3.5 

7.1 

60 

4.6 

9.4 

70 

7.1 

14 . 1 

75 

10.5 

14 . 1 

80 

14 . 1 

14 . 1 

>  81.1 

14.1 

14 . 1 

Denoting  pixel  beginning  and  end  by  subscripts  1  and  2, 
respectively,  the  IFOV  is  then  ^6  -  . 

Solution  of  Equation  6  for  and  yields  an  ECA 
value  of  A*  =  *2  -  .  Referring  to  Figure  2,  the  along- 

line  resolution,  aL,  is  approximately  RA4'/cos7  where  y  is 
the  scan  skew. 

Thus,  the  AVHRR  IFOV  of  1.3  mr  corresponds  to  HRPT 
nadir  resolutions  of  1.08  km  and  1.13  km  for  the  nominal 
satellite  heights  of  833  km  and  870  km,  respectively. 
Along-line  resolution  decreases  progressively  with  increas¬ 
ing  scan  angle  as  illustrated  in  Figure  15  for  a  nominal 
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Figure  15.  HRPT  Along-Line  Resolution 

satellite  height  of  870  km  with  earth  radius  of  6378.14  km. 

APT  processing  reduces  resolution  along  the  scan  line 
by  averaging  HRPT  pixels.  Because  of  IFOV  overlap  with 
consecutive  pixels  (see  Figure  6),  the  IFOV  of  an  APT  pixel 
averaged  from  n  HRPT  pixels  is 

IFOV^pt  =  (n  -  1)STEP  +  IFOV^^p^  (15) 

where  STEP  is  the  mirror  step  size  for  each  AVHRR  word. 
Recall  that  each  step  is  25  timing  counts,  so  that 
STEP  =  0.94439882  mr. 

The  HRPT/APT  reduction  algorithm  in  Table  7  produces  a 
nearly  uniform  resolution  of  4  km  along  the  scan  line  as 
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illustrated  in  Figure  16.  The  AVHRR  IFOV  of  1.3  mr  cor¬ 
responds  to  APT  subpoint  resolutions  of  3.44  km  and  3.59  km 
for  the  nominal  satellite  heights. 

Because  of  the  varying  size  of  HRPT  pixel  images  along 
the  scan  line,  some  sort  of  data  resampling  is  required  for 
HRPT  to  remove  distortion.  Repeating  pixels  in  inverse 
proportion  to  their  resolution  is  one  method,  as  presented 
by  Emery  [3:1177].  Since  AFT  scan  lines  have  nearly  uniform 
resolution,  data  resampling  should  not  be  required  for  APT. 

Subsequent  HRPT  scan  lines  are  nominally  contiguous  at 
the  subpoint  [15:43].  The  line  width  is  given  by  the  for¬ 
ward  progress  of  the  satellite  during  one  scan  line. 

For  a  nominal  satellite  period  of  102  min  and  AVHRR  scan 
rate  of  6  lines/sec,  this  corresponds  to  an  EGA  of 
0.00980  deg  (0.171  mr) .  Using  Equation  13,  the  HRPT  line 
width  is  about  1.09  km. 

Similarly,  the  APT  line  rate  of  2  lines/sec  corresponds 
to  an  APT  line  width  of  about  3.27  km.  It  should  be  noted 
that  each  APT  line  is  actually  the  data  from  one  in  three 
AVHRR  lines.  The  other  two  lines  are  missing  in  APT.  If 
the  APT  lines  are  represented  as  contiguous,  there  is  an 
implicit  repetition  of  one  AVHRR/APT  line  three  times. 

Denoting  line  width  as  along-track  resolution,  the  APT 
and  HRPT  subpoint  resolutions  are  summarized  in  Table  15. 
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Figure  16.  APT  Along-Line  Resolution  [15:17] 


Table  15.  HRPT/APT  Subpoint  Resolution 


Along-Line  fkm) 
HRPT  1.08-1.13 

APT  3.44-3.59 


Alona-Track  (km) 
1.09 
3 . 27 


The  along-track  resolution  is  also  reduced  with  increasing 
scan  angle  due  to  the  increasing  range  from  satellite  to 
image  area  and  the  earth's  curvature. 
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Earth  Location  of  Pixel  Images 

The  earth  scene  viewed  in  each  pixel  represents  a 
square  view  of  a  spheroid  earth,  generally  from  a  skewed 
viewing  angle.  In  general,  the  area  viewed  will  have  an 
irregular  shape  that  is  different  for  each  pixel. 

Here,  the  location  of  a  pixel  will  be  defined  as  the 
earth  location  in  the  center  of  the  pixel  IFOV.  Assuming 
that  the  earth  scene  is  essentially  flat  in  the  pixel  IFOV, 
this  location  also  corresponds  to  the  center  of  the  earth 
area  viewed.  Each  pixel  is  then  associated  with  a  unique 
earth  location.  The  following  chapter  shows  how  to  deter¬ 
mine  the  earth  location  of  each  pixel  and  to  find  the  pixel 
that  contains  a  view  of  a  given  earth  location. 
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Introduction 


Spatial  registration  relies  on  the  fact  that  each  pixel 
is  directly  related  to  an  earth  location  by  the  satellite 
position,  earth  position,  and  sensor  scan  angle.  Knowing 
these  relations  as  functions  of  time  allows  transformation 
from  pixel  coordinates  to  geographic  coordinates  (direct 
referencing)  and  vice  versa  (inverse  referencing). 

The  most  critical  factor  in  spatial  registration  is 
timing.  At  an  AVHRR  sensor  scan  rate  of  six  lines  per 
second,  a  timing  error  of  one  second  in  the  computed  satel¬ 
lite  position  is  equivalent  to  an  error  of  six  AVHRR  lines 
(6.54  km  at  the  satellite  subpoint).  To  account  for  this, 
provision  should  be  made  for  a  time  adjustment  in  the  refer¬ 
ence  time,  tg  (see  Equation  46) . 

To  make  this  time  adjustr  2nt  prior  to  direct  registra¬ 
tion,  the  inverse  registration  algorithm  can  be  used  to 
overlay  a  geographic  outline  on  the  satellite  image.  The 
reference  time  can  then  be  adjusted  to  improve  the  registra¬ 
tion  visually.  Once  a  more  accurate  reference  time  is 
obtained,  the  entire  image  can  be  accurately  registered 
using  the  direct  referencing  algorithm. 
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First,  this  chapter  explores  the  general  relation¬ 
ships.  Second,  various  coordinate  transformations  necessary 
to  relate  different  coordinate  systems  are  presented. 
Finally,  specific  registration  algorithms  are  introduced  for 
both  direct  and  inverse  image  referencing. 

Notation 

In  this  chapter,  the  following  notation  is  observed: 
r  =  vector 

r  =  scalar  magnitude  of  vector  r 
r  =  unit  vector  in  direction  of  r  (r  =  r/r) 
rjj,  ry,  =  scalar  components  of  r  (r  =  xr^  +  yry  +  zr^) 
r  X  y  =  vector  cross  product  of  r  and  y 
r  .  y  =  scalar  dot  product  of  r  and  y 

All  vectors  are  referenced  to  an  inertial  coordinate 
system  (Earth  Centered  Inertial)  with  origin,  0,  at  the 
center  of  the  earth,  and  unit  vectors  defined  as  follows: 

X  =  direction  from  O  toward  the  first  point 
of  Aries  (star  reference) 

y  =•  direction  from  O  in  the  plane  of  the  earth's 
equator  such  that  x  x  2  =  z 
(right-hand  coordinate  system) 

z  =  direction  from  O  toward  the  earth's  north  pole 
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General  Relationships 

In  the  Earth  Centered  Inertial  (ECI)  coordinate  system 
described  above,  the  pixel  earth  location  is  determined  by 
the  satellite  velocity  and  range  vectors  by  the  following 
(see  Figure  17) : 

R  =  r  d  (16) 

where 

R  =  earth  location  vector 
r  =  satellite  position  vector 

d  =  range  vector  from  satellite  to  earth  location 

All  three  vary  in  time.  If  the  satellite  position  r  is 
known  as  a  function  of  time,  then  given  either  R  or  d  at  a 
given  time,  the  other  can  be  computed  with  Equation  16. 
Define  the  scan  vector  §  in  Figure  17  by 

s  =  r  X  y  (17) 

where  y  is  the  satellite  velocity  vector.  Note  that  the 
unit  vector  s  represents  the  scan  direction  (right  to  left) . 

The  off-nadir  scan  angle  5  determines  the  range  direc¬ 
tion  d  from  r  and  s: 

d  =  -rcos4  +  ssin5  (18) 

Recall  that  Equation  3  gives  the  range  distance  d; 

d  =  rcos<S  -  (R  ^  -  r^sin^6)  (19) 

Thus,  combining  Equations  18  and  19  with  the  vector  identity 
d  =  dd  yields  an  expression  for  d  in  terms  of  5: 
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Figure  17.  Satellite/Earth  Geometry  (ECI  Coordinates) 
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d  =  -  r { cosf  [ rcos5  -  (R  ^  ]) 

+  s,{  sinS  [rcos(S  -  (R  ^  -  r  ^sin  ^5) ])  (20) 

Spherical  Earth.  If  a  spherical  earth  is  assumed,  the 
earth  radius  R  is  a  known  constant.  Substituting  Equation 
20  into  Equation  16  then  gives  the  desired  earth  location  R. 

Oblate  Spheroid  Earth.  For  an  oblate  spheroid  earth, 
however,  the  earth  radius  R  is  a  function  of  latitude.  An 
expression  for  R  is  therefore  required  in  terms  of  the 
satellite  radius,  velocity,  and  scan  angle. 

Each  point  (Rj^,  Ry,  R2)  on  an  earth  meridian  ellipse 
satisfies  the  relation 

(Rx"  +  1  (21) 

where  R^^  and  Rp^  are  the  earth  equatorial  and  polar  radii, 
respectively.  Substituting  the  Pythagorean  relation 
Rjj^  +  Ry^  =  R^  -  R2^  into  Equation  21  and  solving  for  R^ 
gives 


where 

«  -  1  (21) 

Note  that  f  =  e^/(l  -  e^  )  where  e  is  the  eccentricity  of 
the  meridian  ellipse. 

From  Equation  16,  R2  =  r2  +  d2  =  Tg  +  ddg  .  Substitu¬ 
tion  into  Equation  22  yields 

-  Urz  +  daz)^  (24) 
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Recall  from  Equation  1  that 

-  2rdcos(S  (25) 

Combining  Equations  24  and  25  then  rearranging  produces 
d^  (1  +  ?d2^  )  -  2d(rcos5  -  ^^2^2)  =  -  r^  -  Cr2^  (26) 

Solving  for  d  by  taking  the  negative  root  (closer  range) , 
gives 

d  =  rcoS(S  -  Cr2d2/(1  +  ?d2^)-  { [Req^  - 

-  ^r2^  +  (rcos5  -  Cr2d2)^]/(1  +  )  )^'^  (27) 

The  range  vector  is  again  derived  by  applying  the  relation 
d  =  dd  to  Equations  18  and  27. 

Satellite  Position 

For  these  registration  methods,  highly  accurate  satel¬ 
lite  positions  are  requ'red.  This  section  will  briefly 
discuss  the  general  orbital  case  and  a  local  approximation 
to  simplify  the  numerical  calculations. 

Elliptical  Orbits.  In  general,  the  satellite  will 
follow  an  elliptical  orbit  with  one  focus  at  the  center  of 
the  earth.  Perturbations  caused  by  the  earth's  equatorial 
bulge,  drag,  and  the  gravitational  attractions  of  the  sun 
and  moon  will  alter  the  orbit  over  time.  To  get  an  accurate 
satellite  position  at  a  time  after  epoch  (time  the  satellite 
position  was  measured) ,  an  elliptical  orbit  model  with 
perturbations  should  be  used.  The  registrations  in  this 
work  used  satellite  position  calculations  produced  by  an 
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elliptical  orbit  program  based  on  the  NORAD  SGP  (Simplified 
General  Perturbations)  model. 

Spherical  Orbit  Approximations.  Over  small  periods  of 
time,  the  satellite  motion  may  be  approximated  by  a  spheri¬ 
cal  orbit.  The  following  discussion  shows  how  approximate 
satellite  position  and  velocity  vectors  may  be  calculated  in 
the  vicinity  of  an  initial  set  of  accurate  satellite  posi¬ 
tion  and  velocity  vectors. 

For  a  satellite  in  a  circular  orbit  at  radius  r,  the 
satellite  circular  velocity,  v^  ,  is  given  by  [2:34] 

V(.  =  (28) 

where 

H  =  earth  gravitational  constant 

(3.98600448  X  10^^  m^  s'^  [20:K6]). 

A  local  approximation  for  the  satellite  angular  velocity,  9, 

is  obtained  by  combining  Equation  28  with  the  relation  v^ 

o 

=  re  : 

8  =  (29) 

For  a  spherical  orbit,  the  angular  speed  is  simply 

8  =  27r/P  (30) 

where  P  ’s  the  orbit  period. 

The  circular  velocity  vector,  y^.  ,  is  given  by 

Yc  =  (-r  X  h)Vc  (31) 

where 

r  =  satellite  radius  vector 

y  =  satellite  velocity  vector 
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h  =  angular  momentum  vector  (r  x  v) 

V(.  =  circular  velocity  from  Equation  28 

Given  the  satellite  position  r(t)  and  velocity  y(t)  at 
time  t,  the  satellite  position  vector  r(t  +  At)  and  veloc¬ 
ity  vector  y(t  +  At)  at  time  t  +  At  can  be  approximated 
using  a  circluar  orbit.  It  can  be  shown  that  they  are  given 
by  [2:72] 

r(t  +  At)  =  r (t)  [cos (^At)  ]  +  v^.  (t)  [sin(5At)/^  ]  (32) 

and 

y^.  (t  +  (St)  =  -r(t){[V(.  (t)/r]sin(^At)  } 

+  Yc  (t) [cos ( At) ]  (33) 

where 

y^  (t)  =  circular  velocity  vector  given  by  Equation  31 
d  =  satellite  angular  velocity  given  by  Equation  29 

The  above  equations  can  be  used  in  the  direct  referencing 
algorithm  to  calculate  satellite  and  subpoint  positions  more 
quickly,  instead  of  using  elliptical  model  calculations  for 
every  pixel. 

Pixel  Time 

The  general  relations  above  show  how  the  earth  location 
vector  is  determined  by  the  satellite  position  and  velocity 
vectors  and  the  sensor  scan  angle  at  a  given  time.  This  is 
sufficient  for  direct  referencing  since  the  pixel  time  is 
known. 
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For  inverse  referencing,  the  earth  location  is  known, 
but  not  the  time  of  the  pixel  which  views  that  location. 

This  time  is  found  using  an  iterative  procedure  by  estimat¬ 
ing  the  time  that  the  scan  vector  will  cross  the  earth 
location  in  question,  assuming  no  earth  rotation.  Updating 
the  earth  position  (in  static  coordinates)  by  accounting  for 
earth  rotation,  and  repeating  the  process  produces  increas¬ 
ingly  better  time  estimates. 

For  the  purpose  of  simplification,  assume  that  the 
earth  is  spherical  and  the  satellite  orbit  circular.  Refer¬ 
ring  to  Figure  18,  let  the  earth  location  and  satellite 
subpoint  geocentric  coordinates  (a,#)  be  denoted  by  the  sub¬ 
scripts  E  and  S',  respectively. 

The  first  step  is  to  find  the  equivalent  time  and 
longitude  of  the  ascending  node  N  (equator  crossing) ,  as¬ 
suming  a  spherical  orbit.  If  the  orbit  is  truly  circular, 
these  will  be  the  actual  values. 

Given  the  geocentric  latitude  fj,  and  longitude  Aj,  for 
subpoint  S'  at  time  t^,  (the  center  of  the  image  is  good), 
Equation  9  can  be  rewritten  to  give  the  satellite  orbit 
angle,  9^,,  from  equator  crossing; 

sintfj,  =  sin*5,/sin(i)  (34) 

where  °9  =  orbit  angular  velocity,  given  by  Equation  3  0 
(spherical  orbit)  or  Equation  29  (spherical  approximation) . 
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Figure  18.  Satellite  Track/Scan  Geometry  (Static) 
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The  time  of  the  ascending  node,  t^  ,  is  then 

(35) 

The  subpoint  static  longitude,  A '5,  is  given  by 

cosCA'j,)  =  cos^g,/cos*5,  (36) 

The  longitude,  A^  ,  of  the  ascending  node  N  can  now  be  com¬ 
puted  as  follows: 

\  =  Ag.  +  A'g,  -  (tg.  *  t^  )  (A  "  A)  (37) 

where 

Q  =  earth  rotation  rate  (0.2506845  deg/min) 

A  =  orbital  precession  rate  (6.844774  x  10’^  deg/min 
for  sun-synchronous  orbits) 

Given  an  earth  location  E  at  geocentric  coordinates 
(♦e  )'  next  step  is  to  guess  a  viewing  time  t^  (use 

the  subpoint  time,  tg,,  or  any  central  time).  The  static 
longitude  of  E,  ,  is  then  given  by 

^'e  =  "  ^e  )  (J^  -  A)  (38) 

Values  for  )3,  j,  and  ^  are  then  computed  using  trigo¬ 
nometric  relations  in  Figure  18: 

COS0  =  COS*gCOSA'g  (39) 

sin(j)  =  sin*g /sin/9  (40) 

sin#  =  sin(  j-i)  sin;0  (41) 

where 

/3  =  great  circle  arc  from  ascending  node  N  to 
earth  location  E 

j  =  angle  between  the  equator  and  the  great  circle 
through  N  and  E 

#  =  earth  center  angle  between  subpoint  S '  and 
earth  location  E 
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Let  S"  be  the  projected  subpoint  location  where  the 
scan  vector  s  crosses  earth  location  E.  The  satellite  orbit 
angle,  9^„,  at  subpoint  S'*  is  given  by  applying  the  cosine 
law  in  triangle  ENS”: 

cos^5„  =  cosyS/cos*  (42) 

The  projected  scan  time,  tj,,,  is  then 

tj..  = 

During  the  time  between  tg,,  and  t^  ,  however,  the  earth 
location  actually  moves  (in  static  coordinates)  to  position 
E"  (see  Figure  18)  .  The  projected  time  tg,,  can  then  be  im¬ 
proved  by  estimating  a  new  scan  crossing  time,  t',  and 
repeating  the  calculations  above.  The  new  estimate  t'  is 
chosen  as 

t'  =  tg  +  a(tg,.  -  tg  )  (44) 

where  a  is  a  convergence  constant  chosen  to  optimize  conver¬ 
gence  efficiency. 

For  an  inclination  over  90  degrees,  assuming  no  earth 
rotation  tends  to  overestimate  the  crossing  time.  A  good 
value  for  a  is  given  by  the  empirical  formula 

a  =  1  +  cos(i)cos*g  (n  -  ^)/°9  (45) 

In  practice,  using  a  value  of  a  =  1  performed  nearly  as 
well . 

This  process  is  repeated  until  the  time  difference  is 
less  than  some  arbitrary  value,  e;  i.e.  |  tg„  "  tg  |  <  e  , 
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This  results  in  an  HRPT  line  accuracy  of  e/T^  lines,  where 
is  the  HRPT  line  period  (1/6  sec) . 

The  process  above  produces  the  best  time  estimate  under 
the  assumption  of  a  circular  orbit  in  the  vicinity  of  sub- 
point  S'.  If  greater  accuracy  is  desired,  the  actual  posi¬ 
tion  of  subpoint  S"  at  time  tj.,  should  be  computed  using  an 
elliptical  orbit  model.  The  iterative  process  of  estimating 
the  scan  time  is  then  repeated. 

A  slight  error  remains  due  to  the  oblate  (flattened) 
shape  of  the  earth.  If  necessary,  this  could  be  corrected 
by  computing  the  pixel  geocentric  coordinates  at  the  esti¬ 
mated  time  using  the  direct  referencing  algorithm  for  an 
oblate  spheroid  earth.  The  time  estimate  could  then  be 
further  improved  through  iteration,  comparing  the  computed 
earth  location  to  the  desired  one. 

Coordinate  Transformations 

The  general  relations  above  reference  Earth  Centered 
Inertial  (ECI)  coordinates.  For  spatial  registration,  pixel 
and  geographic  coordinates  are  desired.  The  following  sec¬ 
tions  show  how  to  transform  between  these  coordinate  sys¬ 
tems  and  ECI  coordinates. 

Pixel  Coordinates.  The  HRPT  data  consists  of  a  series 
of  pixels  referenced  by  integer  line  and  pixel  numbers 
(l,p).  Pixel  numbers  range  from  0  to  2047,  and  line  numbers 
from  0  to  around  5400. 
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Line  and  pixel  coordinates  (l,p)  can  be  considered  as 
continuous  to  simplify  the  mathematical  transformation  to 
and  from  continuous  geographic  coordinates.  The  corner  of 
the  first  pixel  in  the  first  line  has  pixel  coordinates 
(-0,5, -0.5)  with  its  center  at  (0,0).  Pixel  (l,p)  covers 
line  values  (1  -  0.5  to  1  +  0.5)  and  pixel  values  (p  -  0.5 
to  p  +  0.5). 

Define  the  scan  angle  <S  as  negative  at  the  beginning  of 
a  scan  line,  zero  at  nadir,  and  positive  thereafter.  The 
scan  angle  <Sp  at  the  center  of  pixel  (l,p)  is 

=  STEP(p  -  1023.5)  (46) 

where  STEP  is  the  mirror  step  angle  for  each  pixel 
(0.9439882  X  10'^  rad)  . 

The  time  t^  p  corresponding  to  pixel  (l,p)  is  given  by 
ti,p  =  to  +  T,  (1  -  Iq  )  +  Tp  (p  -  Po  )  (47) 

where 

tg  =  reference  time  of  pixel  (Ig  ,Po  ) 

T^  =  time  duration  of  each  line  (51.27778  x  10'^  sec) 

Tp  =  time  duration  of  each  pixel  (0.0250401  x  10'^  sec) 

Conversely,  given  scan  angle  S  and  time  t  for  a  par¬ 
ticular  image  point,  the  line  and  pixel  coordinates  are 
given  by 

p  =  1023.5  +  (5/STEP  (48) 

1  =  Ig  +  [t  -  tg  -Tp  (P-Pg  )  ]/T^  (49) 

Coordinates  1  and  p  will  be  integers  only  if  the  image  point 
happens  coincide  with  the  center  of  pixel  (l,p)  at  time  t. 
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In  general,  it  will  not.  To  find  the  pixel  containing  the 
image  point,  1  and  p  are  rounded  to  the  closest  integer. 
This  assumes  negligible  earth  rotation  (7.292116  x  10'^ 
rad/sec)  during  one  scan  line  (1.215353  x  10'®  rad/line  or 
0.075  km/ line  at  the  equator). 

Valid  pixel  coordinates  are  -0.5  to  2047.5  for  p,  and 
-0.5  to  l^j^  +  0.5  for  1.  Coordinates  outside  this  range 
represent  earth  locations  not  viewed  in  the  data. 

To  register  APT  data,  a  transformation  between  APT  and 
HRPT  pixel  coordinates  is  used  as  summarized  in  Table  16. 
Note  that  the  center  of  APT  pixel  (l,p)  also  has  coord¬ 
inates  (l,p)  and  covers  the  range  (1  +  0.5,  p  +  0.5). 

To  used  sampled  APT  data,  the  pixel  number  must  be 
converted  to  its  data  sample  number.  The  APT  data  rate  is 
4160  words/sec.  If  the  data  is  digitized  at  9600  sam¬ 
ples/sec,  the  pixel  number  should  be  multiplied  by  2.30877 
(9600  samples/4160  words) . 

Inertial/Geocentric  Coordinates.  It  is  convenient  to 
reference  earth  locations  by  latitude  and  longitude  co¬ 
ordinates.  Denote  geocentric  latitude  and  longitude  by  ♦gj. 
and  Agj.,  respectively.  These  can  be  expressed  in  terms  of 
the  components  of  the  ECI  earth  position  vector  R  [10:C- 
13]  : 

tan*g^  =  [R^  /(R^2  +  Ry^  )’^2  j 

tan[Ag^  +  GHA(t)]  =  tan(Ry  /R^)  (51) 
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Table  16.  APT/HRPT  Pixel  Coordinate  Transformation 


Reoion 

HRPT  pixels 

APT 

pixels 

Transformation 

* 

5 

0 

to 

120 

0 

to 

120 

Ph 

= 

Pa  (  1  ) 

4 

121 

to 

213 

121 

to 

182 

Ph 

= 

Pa  (1-5) 

- 

60.5 

3 

214 

to 

379 

183 

to 

265 

Ph 

- 

Pa  ( 

2  ) 

- 

162 

2 

380 

to 

709 

266 

to 

375 

Ph 

= 

Pa  ( 

3  ) 

- 

418 

1 

710 

to 

1337 

376 

to 

532 

Ph 

= 

Pa  ( 

4  ) 

- 

794 

2 

1338 

to 

1667 

533 

to 

642 

Ph 

= 

Pa  ( 

3  ) 

- 

261 

3 

1668 

to 

1833 

643 

to 

725 

Ph 

= 

Pa  ( 

2  ) 

+ 

382 

4 

1834 

to 

1926 

726 

to 

787 

Ph 

= 

P,  (1.5) 

+ 

745 

5 

1927 

to 

2047 

788 

to 

908 

Ph 

= 

Pa  ( 

1  ) 

+ 

1139 

Channel 

A:  1, 

H  ~ 

1a  ( 

3  )  - 

1 

Channel 

B:  1 

H  “ 

1a  ( 

3  ) 

*  =  HRPT  pixel;  p^  -  APT  pixel 

1^  =  HRPT  line;  1^  =  APT  line 

Note;  Because  of  the  overlapping  of  APT  pixels  in 
region  4,  the  pixel  centers  are  offset  from  the  IFOV 
center  by  0.25  STEP  in  region  4. 


v.here  GHA(t)  the  Greenwich  hour  angle  at  time  t.  This  is 
computed  in  reference  to  the  beginning  of  the  year  by  [10:C- 
12;  2:99-100] 

GHA(t)  =  Gq  +  G,  t^  +  ntf,  (52) 
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where 


Gg  =  GHA  at  the  beginning  of  the  year  of  time  t  (deg) 
G^  =  daily  increase  in  GHA  (0.985647366  deg/day) 
t^  =  day  of  the  year  of  time  t 
n  =  rate  of  earth  rotation  (360.9856476  deg/day) 
tf,  =  fractional  part  of  a  day  in  time  t 

Note  that  GHA  is  measured  from  the  Greenwich  meridian  west¬ 
ward  to  the  first  point  of  Aries.  Values  of  Gg  can  be 
found  in  the  Astronomical  Almanac  [20]  for  the  current  year 
or  calculated  in  reference  to  epoch. 

Geocentric/Geodetic  Coordinates.  Map  coordinates  are 
usually  geodetic  (also  called  geographic) .  As  illustrated 
in  Figure  19,  geodetic  latitude  is  based  on  the  local  hori¬ 
zon  of  an  oblate  earth,  while  geocentric  latitude  is  based 
on  the  center  of  the  earth.  Geodetic  longitude  is  the  same 
as  geocentric  longitude. 

Denoting  geodetic  latitude  and  longitude  by  and 
respectively,  geodetic  coordinates  are  related  to  geocen¬ 
tric  coordinates  by  [3:1177;  10:C-10] 

=  (R«/^po)' tan*g,  (53) 

V  =  V 

where 

=  Ecjuatorial  radius  of  the  earth  (6378.144  km) 

Rp^  =  Polar  radius  of  the  earth  (6356.759  km) 
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Horizon 


S 


Legend ; 


E:  Earth  station 
0:  Earth  center 

♦ggt  Geocentric  latitude 
Geodetic  latitude 


Figure  19.  Geodetic  vs.  Geocentric  Latitude 


Direct  Referencing 

In  direct  referencing,  line  and  pixel  coordinates  (l,p) 
are  converted  to  geodetic  coordinates  for  a  given 
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pixel.  Once  the  coordinates  of  each  pixel  are  obtained,  the 
data  can  be  displayed  using  any  projection  method  desired, 
resampling  as  required. 

Assuming  an  oblate  spheroid  earth,  a  direct  referen¬ 
cing  algorithm  is  presented  below. 

STEP  0:  Perfomn  initial  setup. 

0.1  Calculate  reference  time  tg  for  reference  pixel 
(Ig  ,Pg  )  using  mark  time  t  at  pixel  coordinates 
(l,p)  in  Equation  47  (use  HRPT  pixel  coordinates) 

0.2  Calculate  satellite  position  r  and  velocity  y  at 
reference  time  tg  using  a  general  elliptical 
orbit  model 

0.3  Calculate  subpoint  geocentric  latitude  and 
longitude  Aj.  using  Equations  50  -  52 

0.4  If  using  spherical  orbit  approximation,  calculate 
ascending  node  longitude  and  time  t|j  using 
Equations  34  -  38. 

For  each  pixel  with  coordinates  (l,p),  do  the  following: 

STEP  1:  Calculate  the  pixel  time  and  scan  angle 

1.1  If  APT  pixel  coordinates,  transform  to  HRPT  pixel 
coordinates  using  Table  16 

1.2  Calculate  the  pixel  scan  angle  6  with  Equation 
46 

1.3  Calculate  the  pixel  time  t^ ^  with  Equation  47 

STEP  2:  Compute  satellite  position  vector  r  and  velocity 
vector  y  at  time  t^  p 

2.1  It  using  a  general  orbit,  compute  r  and  y  using  an 
elliptical  orbit  model 

2.2  If  using  a  local  spherical  orbit  approximation, 
compute  £  and  y  using  Equations  32  and  33 

STEP  3 :  Compute  or  update  scan  direction  s  as  required 
using  Equation  17 
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STEP  4: 

Compute  range  vector  d  from  satellite  to  image 
point 

4 . 1 

Calculate  range  direction  d  using  Equation  18 

4.2 

Calculate  range  distance  d  using  Equation  19  for  a 
spherical  earth,  or  using  Equation  27  for  an 
oblate  spheroid  earth. 

4 . 3 

Calculate  range  vector  d  with  the  identity  d  = 

dd 

STEP  5: 

Compute  earth  position  vector  R  using  Equation 

16 

STEP  6: 

Calculate  geocentric  latitude  and  longitude 

using  Equations  50  and  51 

STEP  7: 

Convert  to  geodetic  latitude  ♦gj  and  longitude 
using  Equations  53  and  54 

^gd 

Inverse  Referencing 

For  inverse  referencing,  the  object  is  to  find  the 
pixel  that  views  a  given  earth  location.  This  involves  an 
iterative  procedure  of  estimating  the  time  of  viewing,  then 
making  corrections  for  earth  rotation  to  get  progressively 
better  time  estimates. 

Once  the  viewing  time  is  known  (with  sufficient  accur¬ 
acy)  ,  the  pixel  coordinates  in  the  image  can  be  calculated. 
This  can  be  used  to  overlay  geographic  outlines  or  grids  on 
the  satellite  images. 

The  inverse  referencing  algorithm  that  follows  applies 
to  both  elliptical  orbits  and  spherical  orbit  approxima¬ 
tions  . 
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STEP  0: 


Perforin  initial  setup. 


For  a 

STEP 

STEP 


STEP 

STEP 


.1  Calculate  reference  time  tg  for  reference  pixel 
(1q  ,pg  )  using  mark  time  t  at  pixel  coordinates 
(l,p)  in  Equation  47  (use  HRPT  pixel  coordinates) 

.2  If  using  a  spherical  orbit  approximation,  calcu¬ 
late  satellite  position  r  and  velocity  v  at  refer¬ 
ence  time  tg  using  a  general  elliptical  orbit 
model 

given  earth  location  E  with  geodetic  coordinates 
,)  ,  do  the  following: 


1:  Convert  from  geodetic  coordinates  to  geocentric 

latitude  and  longitude  A  for  earth  location 
E,  using  Equations  53  and  54 

2:  Compute  the  equivalent  time  and  longitude  of  the 

ascending  node  (equator  crossing) 

.1  Calculate  a  reference  subpoint  position  (ECI 

coordinates)  at  a  reference  time  (say  tg  )  using  a 
general  elliptical  orbit  model 

. 2  Calculate  the  geocentric  latitude  *5,  and  lon¬ 
gitude  Aj,  for  subpoint  S'  with  Equations  50  and 
51 

.3  Use  Equations  34-37  to  calculate  the  time  t^  and 
longitude  A^  of  ascending  node  N 

3:  Guess  a  viewing  time  tg  (tg  is  good) 

4:  Estimate  the  scanning  time  of  E,  assuming  no  earth 

ro-*  Ation 

.1  Compute  the  static  longitude  of  E,  A'g  ,  using 
Equation  38 

.  2  Compute  related  geometric  measures  ^ ,  j  ,  and 
using  Equations  39,  40,  and  41,  respectively 

.3  Using  Equations  42  and  43,  calculate  the  estimated 
scanning  time  t^„  when  the  satellite  reaches  sub- 
point  S” 
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STEP  5: 


5.1 


5.2 
STEP  6: 

6.1 

6.2 

STEP  7: 
7 . 1 


7.2 


7.3 
STEP  8; 
STEP  9: 


9.1 


9.2 


Make  an  improved  scan  time  guess,  t',  given  by 
Equation  44 

For  better  convergence,  use  convergence  factor  a 
computed  using  Equation  45 

Using  a  =  1  is  easier,  and  nearly  as  good 


Check  the  difference  between  tj,,  and  tg 


If 

^  t  —  t  ^ 
1  '“S"  E  1 

<  €,  go  to  STEP  7 

If 

*  t  “  t  ^ 

1  ^S"  E  1 

>  €,  return  to  STEP  4..  substituting 

t ' 

for  t 

Consider  assumptions  and  accuracy  desired 

If  local  spherical  orbit  assumption  is  unsatisfac¬ 
tory,  return  to  STEP  2,  using  an  elliptical  model 
to  calculate  the  projected  position  of  subpoint  S' 

If  spherical  earth  approximation  is  unacceptable, 
compute  earth  location  using  direct  referencing 
algorithm  with  an  oblate  spheroid  earth,  and 
compare  to  actual  earth  location  to  improve  time 
estimate 

Otherwise,  proceed  to  STEP  8 

Compute  scan  angle  from  ^  using  Equation  7 

Compute  pixel  coordinates  (l,p) 

Calculate  HRPT  line  number  1  and  pixel  number  p 
using  Equations  48  and  49,  rounding  to  the  closest 
integer 

If  using  APT  data,  convert  from  HRPT  to  APT  line 
and  pixel  numbers  using  Table  16 
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V.  Conclusions /Recommendations 


Conclusions 

All  the  information  necessary  for  the  spatial  registrar 
tion  of  TIROS-N  weather  satellite  images  has  been  presented. 
The  background  theory  has  been  explained  and  developed  to 
meet  registration  requirements  for  both  HRPT  and  APT  data. 

Specific  registration  algorithms  were  presented  to 
convert  pixel  coordinates  to  geographic  coordinates  (direct 
referencing)  and  from  geographic  coordinates  to  pixel  co¬ 
ordinates  (inverse  referencing) ,  These  algorithms  are 
applicable  for  elliptical  orbits,  spherical  orbits,  or  local 
spherical  orbit  approximations.  The  direct  referencing 
algorithm  allows  the  use  of  an  oblate  spheroid  earth  model, 
or  a  spherical  earth  model.  The  inverse  referencing  al¬ 
gorithm  assumes  a  locally  spherical  earth,  with  guidelines 
for  implementing  an  oblate  spheroid  earth  model. 

Preliminary  testing  using  a  personal  computer  verified 
the  basic  operation  of  the  inverse  referencing  algorithm.  A 
computer  program  applying  the  inverse  referencing  algorithm 
to  overlay  a  map  outline  on  a  satellite  image  transmitted 
via  APT  is  being  written,  though  not  completed  as  of  this 
writing. 
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Recommendations 


Having  clearly  set  forth  the  means  and  use  of  spatial 
registration  of  TIROS-N  weather  satellite  images,  many 
useful  applications  and  further  improvements  are  possible. 
These  include  evaluating  and  improving  the  algorithm  perfor¬ 
mance,  calibration  of  the  sensor  data,  and  using  the  regist¬ 
ered  images  for  environmental  (weather)  analysis.  Specific 
recommendations  are: 


1.  Evaluate  and  improve  registration  algorithms 

a.  Perform  a  complete  validation  and  verification  of 
the  registration  algorithms 

b.  Evaluate  algorithm  performance  using  testing  and 
sensitivity  analysis 

2 .  Develop  associated  computer  programs  to  improve  the  use 

of  registered  images 

a.  Analyze  analog-to-digital  conversion  to  get  the 
best  estimate  of  the  original  AVHRR  words 

b.  Develop  a  computer  program  to  perform  calibration 
of  the  sensor  data 

c.  Study  various  map  projection  techniques  to  display 
the  satellite  images 

3.  Use  image  registration  for  various  applications  using 

TIROS-N  weather  satellite  data 

a.  Design  a  geographic  grid  cell  database  for  satel¬ 
lite  data 

b.  Perform  time  series  analysis  of  registered  satel¬ 
lite  images 

c.  Use  satellite  data  for  image  identification  and 
analysis  (such  as  environmental  and  agricultural 
monitoring) 
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